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ABSTRACT 


This  final  report  is  a  result  of  a  study  performed  for  RADC  under 
Contract  Number  F30602-69C -0042 .  The  purpose  of  the  study  was  to  fit  one 
or  more  prior  distributions  to  8  =  MTBF  ■=  Mean  Time  Between  Failure. 

In  particular,  the  objectives  were  three: 

i)  Establish  criteria  for  data  that  would  be  suitable 
for  fitting  prior  distributions  to  0  »  MTBF. 

ii)  Develop  methods  of  fitting  and  fit  one  or  more 
prior  distributions. 

iii)  Perform  robustness  analysis  of  fitted  prior 
distributions. 

It  was  discovered  that  if  the  number  of  identical  equipments  and  number 
of  failures  observed  per  equipment  are  relatively  small,  special  methods  of 
fitting  are  required.  For  the  data  used  in  this  study,  the  inverted  gamma 
distribution  turned  out  to  be  a  good  prior  fit. 
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EVALUATION 


A  thorough  review  of  the  available  literature  on  Bayesian  reliability 
(as  listed  In  Section  11  of  this  report)  will  reveal  that  very  little. 

If  any,  effort  has  been  directed  towards  developing  methods  for  fitting 
prior  distributions  to  empirical  data.  Usually,  priors  are  determined 
through  engineering  Judgment,  personal  experience,  and  wide-scale 
assumptions.  Such  means  of  developing  priors  are  unsatisfactory  if 
Bayes'  methods  are  to  be  used  in  equipment  reliability  demonstration. 

The  objectives  of  this  Phase  I  study  were  to  (a)  establish  data 
criteria  for  use  in  fitting  priors,  (b)  develop  methods  of  fitting 
and  actually  fit  one  or  more  priors,  and  (c)  perform  robustness 
analysis  of  fitted  priors.  Tbs  results  of  Phase  I  show,  for  the  first 
time,  how  priors  can  be  fitted  to  empirical  data,  the  amount  and  type 
of  data  required,  and  the  effect  on  the  posterior  distribution  of  varying 
parameters  of  the  prior  distribution. 

Each  of  the  objectives  of  Phase  I  was  successfully  completed.  Section  5 
shows,  for  three  foams  of  operational  data,  the  minimum  values  of  "n" 

(number  of  equipments)  and  "k"  (number  of  failures  per  equipment)  required 
for  fitting  a  valid  prior  distribution  to  MEHP.  The  results  of  this 
section  provide  a  firm  foundation  on  which  future  data  collection 
programs  for  fitting  priora  can  be  baaed.  Section  4  presents  methods  of 
fitting  prior  distributions  and  Includes  seven  Inverted  0una  priors 
fitted  to  empirical  data  collected  on  seven  different  types  of  equipment. 

These  results  confirm  the  practicality  of  the  Bayes'  method  in  reliability  * 

demonstration,  and  provide  Justification  for  the  use  of  the  inverted 

gamsa  as  tbe  prior  distribution  on  equipment  KTBF.  Section  7  contains 

the  results  of  the  robustness  analysis  performed  to  investigate  the 

effects  of  errors  in  estimating  tbe  scale  and  shape  parameters  of  the 

inverted  gaema  prior  on  tbe  posterior  inverted  geaaa  distribution.  In 

addition,  as  previously  stated.  Section  11  Includes  a  comprehensive, 

7 B- entry  bibliography  of  sources  on  Bayes '  reliability. 

Tbe  methods  developed  in  this  study  were  based  on  "equipment"  level  data* 
thus  it  was  assumed  that  tbe  conditional  distribution  of  tims-to- failure 
was  assumed  to  be  exponential.  However,  it  should  be  noted  that  the 
methods  given  here  are  generally  applicable  whatever  the  form  of  the 
conditional  distribution. 

Phase  I  results  will  serve  as  important  Inputs  to  the  Phase  Q  study 

which  lc  scheduled  to  start  In  February  iy70.  The  methods  developed  and 

data  criteria  established  In  Phase  I  will  be  used  to  fit  additional 

priors  on  other  equipment  types  in  Phase  II.  In  addition,  this  next  < 
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phase  vill  include  the  investigation  of  methods  of  combining  priors  from 
similar,  hut  not  Identical,  equipment,  phase  II  objectives  also  Include 
the  establishment  of  plans  for  developing  and  implementing  Bayesian 
reliability  demonstration  tests,  as  a  prelude  to  the  projected  FI -71 
Phase  HI  study.  _ 
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0.0  SUMK4BI 


Tha  knowledge  of  the  prior  distribution  is  central  to  the  efficient  use 
of  Bayes  methods  in  reliability  estimation  and  demonstration.  This  study 
presents  methods  and  examples  of  fitting  prior  distributions  when  the  data 
available  is: 

i)  n  sample  MTBF's  on  n  identical  equipments. 

ii)  n  observations  on  the  number  of  failures  occurring  in  a  fixed  time 
T  for  n  identical  equipments. 

This  data  must  be  used  since  samples  from  the  prior  distribution  itself, 
i.e.,  the  true  MTBF,  are  unobtainable .  Two  cases  were  considered:  the  family 
of  the  prior  distribution  being  specified  and  unspecified.  In  case  the 
analyst  is  unwilling  to  specify  the  family  the  methods  are  not  wholly  satis¬ 
factory;  the  most  suitable  being  a  method  due  to  Krutchkoff  and  Rutherford 
which  requires  the  prior  distribution  belong  to  the  Pearson  class.  When  the 
family  is  specified  (e.g, ,  inverted  gamma)  methods  are  presented  which  use 
the  known  marginal  distribution  to  fit  the  prior  distribution.  Some  diffi¬ 
culty  is  experienced  when  the  sample  sizes  (K^)  for  the  n  sample  MTBF's  are 
not  all  identical  and  small.  1>.  :  is  case  a  mixed  model  for  the  marginal 
distribution  is  used.  If  the  are  all  large  (roughly  >30)  even  though 
not  identical  the  marginal  distribution  approaches  the  prior  distribution 
rapidly.  If  n  and  jx^|  are  both  large  the  data  is  fitted  directly  to  the 
prior  distribution. 

The  results  are  given  in  the  following  table : 


General  Methods  of  Fitting  Prior  Distributions 


Family  Unspecified 

(all  values  of  n  and  j  K^J  ) . 

Estimate  first  four  moments 
of  prior  distribution  from 
data  and  use  Pearson  class. 

Family  Specified 

jKj  large,  n  2 

large  enough  to  use  X  te3t. 

|KjJ  small,  n 
large  enough  to  use  mixed 
model  test. 

Treat  data  as  though  it  came 
from  prior  distribution  and 
use  classical  methods. 

Fit  data  to  known  marginal 
using  classical  methods. 

1 


mffv yy 


1.0  INTRODUCTION 


1.1  OBJECTIVE  OF  THE  STUDY 


There  are  a  number  of  measures  of  reliability.  The  most  important  are 

i)  Mean  Time  Between  Failure  (MTBF) . 

ii)  Failure  Rate. 


iii)  Probability  of  Survival  for  a  fixed  time  T. 


iv) 


Time,  say  x  ,  for  which  the  probability  of  survival  is  p.  That  is, 
p  t/h 

x  is  the  (l-p)  quantile  of  the  time  to  failure  distribution. 


It  is  customary  in  equipment/ system  development  to  place  certain  specifica¬ 
tions  on  one  or  more  of  the  above  measures.  Demonstration  tests  (in  statisti 
cal  language  hypothesis  tests)  are  a  fundamental  tool  in  verifying  that  these 
specifications  have  been  met.  Unfortunately,  the  demonstration  tests  sure 
conducted  in  an  environment  of 


•  High  reliability  requirements. 

•  Limited  funds. 


•  Short  time  available. 

Generally,  the  precision  of  the  demonstration  test  increases  as  the  number  of 
failures  observed  increases.  But  high  reliability  requirements  mean  long 
times  to  observe  failures.  On  the  other  hand,  low  producer  and  consumer 
risks  are  desired  and  this  also  means  long  test  times.  This  creates  a  cost/ 
time  problem  tfhich  is  apparently  unsolvable  by  classical  methods.  In  fact, 
Bayes  methods  apparently  hold  the  most  promise  in  solving  the  cost/time 
problem  of  demonstration  tests. 

In  many  sampling  situations,  Bayen  methods  are  not  applicable  because 
the  unknown  parameter  (here,  MTBF,  failure  rate,  etc.)  cannot  be  considered 
a  random  variable,  which  is  a  "must"  to  use  Bayes  methods.  However,  it  is 
clear  that  in  reliability  the  parameter  can  often,  if  not  always,  be  con¬ 
sidered  a  random  variable.  Consider  a  computer  manufactured  by  a  particular 
Company  to  a  particular  design;  each  successive  computer  differing  in  serial 
number  and  parts.  Because  of  these  part  (not  part  type)  differences  and 
other  differences,  each  computer  will  possess  a  different  true  (but  unknown) 
MTBF  **  8  and  hence,  MTBF  *  8  may  be  considered  a  random  variable.  The  essence 
of  Bayes  methods  is  that  a  probability  distribution  is  assumed  to  exist  on  the 
parameter  (here,  measure  of  reliability)  in  question.  This  probability 
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distribution  is  called  the  prior  distribution.* 

The  Measure  of  reliability  considered  throughout  this  study  is  i)  above, 
mutely,  the  very  important  KTBF  -  8.  Not  only  is  MTBF  Most  cossaonly  used  as 
a  Measure  of  reliability,  it  has  the  advantage  that  it  is  a  parameter  in  the 
tiMe  to  failure  distribution,  which  results  in  variables  type  data  rather 
than  ill)  and  iv)  above,  which  result  in  attributes  data.  Thus,  MTBF  permits 
a  parametric  rather  than  a  nonparametric  approach. 

In  susHMry  then,  to  use  Bayes  Methods  a  prior  distribution  is  needed. 

The  basic  objective  of  this  study  is  to  fit  one  or  more  prior  distributions 
to  KTBF  ■  8  for  ground  electronic  equipments.  In  doing  this,  there  are  many 
related  questions  so  that  the  study  objectives  are  mare  particularly  stated 

as: 

i)  Establish  criteria  for  suitable  data  for  fitting  prior  distribu¬ 
tions  to  8r 

il)  Fit  cme  or  more  prior  distributions. 

Hi)  Perform  robustness  analysis  of  fitted  prior  distributions. 

All  these  objectives  were  accomplished  and  are  discussed  in  the  appropriate 
sections. 

The  next  section  gives  More  detail  and  insight  into  the  use  of  Bayes 
Msthods  in  reliability. 

1.2  INTRODUCnClf  TO  BAYES  METHQD6  IN  RELIABILITY 

A  classical  lower  one  sided  confidence  Interval  for  8  (KTBF)  consists  of 

a  statement  that 

P(§  *  8#)  «  1-Y,  Y  small  and  is  prepared  after  the  data  is  available. 
The  probability  can  only  he  interpreted  as  a  confidence.  A  classical  (e.g., 
WL-OTD-781B)  reliability  demonstration  test  consists  of  preselecting  (i.e. , 
before  the  data  is  gathered) 

8*  m  0^  (minimum  acceptable  KTBF)  and  8  *  consumer's  risk.  In  carter  to 

determine  a  unique  test  it  is  also  required  to  "add"  a  8q  (specified  KTBF) 
and  «  »  producer's  risk.  It  I.j  interesting  and  informative  to  dwell  on  the 
"names*  of  8-  and  0q,  i.e.,  minimum  acceptable  and  specified  respectively. 

Okm  might  think,  in  view  of  the  "nasae"  of  8}  that  8q  would  be  called  maximum 
acceptable  or  simply  acceptable  8.  It'  this  were  so,  things  wculd  be  very 
misleading.  That  is,  8q,  8^  are  well  named  nor.  The  name  far  depicts 
what  it  really  la:  a  prespecified  lower  confidence  limit.  The  name  for  8q 
depicts  what  it  really  la:  an  "add  on"  to  obtain  a  unique  test.  Briefly,  it 
is  comaomara  who  daamnd  tests.  They  specify  (8^,0).  Once  having  done  this 
the  producer  selects  tbs  (unique)  test  suitable  to  him. 

In  the  morn  mathematical  literature,  this  prior  distribution  is  sometimes 
called  tbs  Compounding  or  Mixing  distribution. 


k 


i 


Using  Bayes  aethods  It  is  possible  to  prepare  a  probability  (as  against 
confidence)  stateasnt  to  tbe  effect 


P(8  a  e»  (observed  data)  ■  1  -  y 


(1.2.1) 


A  unique  Bayes  deaonstratian  test  can  be  achieved  with  various  "add  ons "  but 
that  is  not  the  subject  of  this  report.  The  iapartant  point  is  the  condition¬ 
ing  rand on  variable  "observed  data"  in  (1.2.1)  above.  It  is  laaaterlal  vhat 
it  is  as  long  as  it  is  a  sufficient  statistic  for  0.  Soaetiaes  the  vector 
of  failure  tines  (x^,  ,  x^)  aay  be  used  and  soaetiaes  the  observed  KTBP 


0  nay  be  used.  In  any  event  to  prepare  (1.2.1),  i.e.,  a  Bayes  confidence 
Interval, 


P(6  2  (observed  data  ■  0)  ■  J  g(0|0)  d0, 

tbe  posterior  distribution  of  0  is  required.  Here  ve  have  illustrated  it 
and  win  continue  to  illustrate  with  the  conditioning  r&ndon  variable  6. 


Clearly, 


g(elf)  -  f(e '»)«(<»  ;  f(S)  t  o. 
/(•) 


(1.2.2) 


Thus,  to  obtain  a  Bayes  confidence  interval  one  aust  have 


i )  The  prior  distribution  g(0). 


11)  Tbe  sampling  distribution  of  tbe  conditioning  statistic  given  0, 


i.e.,  f(0!0). 


Hi)  The  aarglnal  distribution  of  the  statistic  f(t). 


Each  of  these  distributions,  particularly  1)  and  ill)  above,  play  a  vital 
role  in  this  report.  The  object  is  to  fit  g(|).  Tbe  data  gathered  for  this 
study  Is  at  the  equipment  level  and  it  is  assuaed  throughout  that  the  con¬ 
ditional  distribution  of  tlass  to  failure  given  0,  i.e.,  f(xli),  is  exponen¬ 
tial  so  that  f(0l0)  is  os,— .  It  is  then  relatively  staple  to  find  f(f) 
once  the  prior  distribution  is  specified  since 


— 

f($)  -  /  f(4 !•)§(•)  di 


(1-2.3) 


In  eoae  cases  the  observed  randoa  variable  is  the  of  failures  (say  x) 

occurring  in  a  fixed  tlas  period  T  so  that  the  f(x|0)  is  Poisson  and  the 
aarglnal  distribution  of  x  is 


(1.2.4) 


f(x)  -  /  f(x|8)g(0)  d0 
o 

Thu*,  of  the  four  (4)  distributions  of  (1.2.2),  only  two  (2)  are  unknown  in 
the  sense  they  wst  be  estimated  since  the  aargln&l  distribution  in  the 
denominator  is  obtainable  through  g(0)  and  f(dl0). 

1.3  THE  GEHTRAL  PROBLEM 

The  central  objective  of  this  study  is  to  fit  a  prior  distribution  to 
sene  reliability  data;  the  central  problem  is  that  for  a  given  piece  of 
equipment  the  true  0  is  unknown  and  remains  unknown  unless  an  Infinite  number 
of  failures  are  observed.  Thus,  in  fitting  g(0),  one  does  not  even  have 
random  samples  from  g(0).  To  obtain  a  random  sample  from  g(0)  would  be  to 
know  the  MFHF  of  an  equipment  exactly  which  is  Impossible.  There  are  then 
two  sample  sizes  of  concern.  First,  the  number  of  "identical”  equipments 
sampled  which  we  call,  hereafter,  n.  Second,  there  is  the  number  of  failure 
times  available  on  each  equipment  which,  for  the  itb  equipment,  we  call  K* 

and  far  all  n  equipments  we  denote  the  set  of  K.  as  •  The  situation 

is  depicted  below  in  Figure  1.3*  The  squares  represent  identical  equipments. 
In  this  case  x  represents  lifetime. 


Figure  1.3  Variables  Involved  in 

Fitting  a  Prior  Distribution 


At  first  thought,  one  Bight  be  tempted  to  plot  the  8 'a  (somehow  weighted 
far  unequal  K's)  In  a  histogram  and  by  usual  methods  fit  some  prior  distribu¬ 
tion.  To  see  that  this  Is  not  correct,  consider  first  the  case  of  Identical 
K's.  If  the  s' s  were  plotted  in  schist ograa  and  n  taken  very  large  the  histo¬ 
gram  would  look  very  much  like  f^.(0) .  We  are  essentially  "summing"  over  6 

In  the  Joint  density  h^(9,|).  That  Is,  each  sample  Involves  a  pair 
($^,8^)  !■!,  . ..,  n  but  the  0^  component  is  unobservable.  Thus,  each  0  may 
be  considered  a  random  sample  from  f^(8) .  In  the  case  of  unequal  K's  (but  not 

all  necessarily  different)  things  get  ^essy  quickly  but  a  mixed  population 
model  appears  descriptive,  i.e.,  each  8  is  regarded  as  a  random  sample  from 

(9)  “  P-^e)  +  +\  *2^  (1.3.1) 

number  of  Identical  K, ' s 

Where  t  is  the  number  of  distinct  K,  *s  and  p,  ■  _  * 

11  n 

also  Spj  -  1.  The  numerator  in  P1  is  taken  after  renumbering  tc  obtain  different  ^*s. 
This  model  will  be  discussed  In  more  detail  in  Section  4.0. 

The  important  point  of  this  whole  discussion  Is: 

Whatever  methods  are  developed  for  fitting  the  prior  distribution, 
the  data  to  be  used  will  be  data  which  is  taken  from  the  marginal 
distribution  not  from  the  prior  distribution 

The  next  section  gives  definitions  which  will  be  used  throughout  this 
report. 

1.4  amiiiTioifS  akd  assuMPnofs 

Symbol _  Definition 

n  The  number  of  identical  eqjjdfswstt  us&ad  to  fit  a 

prior  distribution. 

Kj  The  number  of  failures  available  on  the  1th  equipment. 

1*1 ,  . . . j  n . 


0  The  unknown  KTHP  of  an  equipment. 


X  The  rondoe  variable  denoting  either  time  to  failure  cc 

washer  of  failures  la  a  fixed  time.  The  context  will 
make  clear  which  is  being  discussed. 

x  A  particular  value  of  X. 


Symbol _ Definition _ 

th 

The  q  quantile  of  the  Cui -Square  distribution  with 
u  degrees  of  freedom. 

The  prior  distribution  on  the  random  variable  t). 

The  posterior  (conditional)  distribution  of  9  having 
observed 

The  conditional  distribution  of  $  for  fixed  9. 

The  marginal  (or  unconditional)  distribution  of  § 
vhere  8  ic  based  on  K  failures. 

The  marginal  (or  unconditional)  distribution  of  time 
to  failure  or  number  of  failures  per  fixed  tame  T. 

The  context  will  make  it  clear  which  is  being  discussed. 

The  conditional  distribution  of  time  to  failure  (or 
failures  per  T)  for  fixed  8. 


Coe  of  the  important  assumptions  of  this  report  is  that  the  conditional 
distribution  of  time  to  failure  for  fixed  8  ic  exponential: 

f(xie)  »  i/&  « ~x/e  x,e  >  0>  (i.4.i) 

Eence ,  the  number  of  failures  occurring  in  fixed  time  T  is  Poisson: 

f(XJe)-S — T,e>o  (1.4.2) 

X  • 

x  =  0,1,  ... 

Tb*  ultimate  aim  of  fitting  prior  distributions  in  reliability  is  Bayec 
demonstration  testa.  8ince  demon strati  on  tests  sue  most  generally  applied 
at  equipment/ syatea  levela  rather  than  part  levels  the  study  was  restricted 
to  this  level  of  data.  Far  this  reason,  the  exponential  assumption  was  felt 
to  be  valid  since  it  has  been  experienced  repeatedly  and  certain  limit 
theorems  (e.g. ,  the  Dr-snick  theorem)  indicate  an  exponential  distribution  is 
to  be  expected  for  complex  equipment.  However,  the  methods  developed  here 
are  generally  applicable  whatever  tbs  form  of  the~"coDditlooal  distribution. 
What  might  happen  though  fear  other  cooditionaidlstributioosio  that  the 
arithmetic  will  rapidly  become  intractable. 


fi 


g(b) 

g(»l$) 

f(£le) 

fR(£) 

f(x) 

f(x|e) 


SECTION  2.0  LATA  COLIECTION 


2.1  TYKES  OF  DATA 

Failure  data  can  occur  in  several  forma.  It  is  divided  Into  tvo  main 
classes:  attributes  and  variables  data.  The  attributes  situation  occurs 
when  an  equipment  is  operated  and  its  survival  or  non  survival  far  tine  T 
(usually  mission  tine)  is  observed.  Due  to  the  nature  of  this  type  of  data, 
attributes  data  was  considered  unusable  far  this  study,  therefore,  no  data  of 
this  type  vas  collected. 

The  variables  data  situation  occurs  when  the  actual  failure  times  are 
available.  These  tines  occur  by  agreeing  to  stop  testing  either  after  a 
fixed  number  of  failures  have  occurred  or  after  a  fixed  tine  has  elapsed. 

The  former  case  Is  called  a  censored  test  and  the  latter  is  called  a 
truncated  test.  Often  the  failure  tines  themselves  axe  not  available  but 
the  pair  (e  »  observed  HTBF,  K  -  number  of  failures)  Is  available. 


2.2  DATA  COLLECTION  PIAN 

The  search  for  potential,  data  sources  vas  initially  limited  to  Bughes 
Aircraft  ttnd  Government  sources.  A  complete  and  comprehensive  search  for 
high  qualit.  reliability  data  on  ground  electronic  equipment  vas  conducted 
in  all  of  Hughes  Aircraft  Divisions  including  the  Quality,  Reliability, 
Effectiveness,  Field  Service  Organisations  and  the  major  Government  data 
centers  maintained  by  the  U.  S.  Government  Agencies.  Much  data  vas  available 
but  due  to  the  high  quality  restriction:  equipment , level  vitb  n  fbirly  large, 
the  amount  of  data  usetul  far  the  establishment  of  a  priori  distributions  vas 
greatly  reduced.  The  data  location  effort  in  the  later  phase  of  the  search 
Included  contacts  in  private  Industry.  IXie  to  the  nature  of  the  data  (KTBF) 
requested  and  proprietary  rights  of  the  contractor  ve  could  obtain  no  data 
from  industry. 

2.3  DATA  COLLECTED 

2-3  1  HUGHES  INTERNAL  SOURCES 

The  primary  source  of  high  quality  data  for  establishing  a  priori  dis¬ 
tributions  vas  Bughes  Aircraft  Company.  The  Bughes  Aerospace  Division  has 
performed  AGREE  tests  on  various  Bughes  built  electronic  equipments.  The 
Hughes  Ground  Systems  Group  has  collected  data  on  Its  systems  and  equipment 
far  many  years.  This  data  includes  both  laboratory  and  field  perfarmt&ce 
data. 


Table  2.3. 1.1  lists  tbs  data  located  and  acquired  vlthiu  Bughes  Aircraft 
for  utilisation  in  this  study. 
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TABLE  2. 3-1*1  HUGHES  AIRCRAFT  COMPANY  ACQUIRED  DATA 


Assigned 
Data  Set 
Number 

Equipment 

n 

EKi 

Type  of  Tests 

9. 

I ram  Computer 

28 

44 

Lab  (AGREE) 

10. 

Multi -mode  Storage  Tube 
Indicator  (Display  Tube) 

20 

24 

Lab  (AGREE) 

11. 

Infrared  Subsystem  in 

Fire  Control  Systems 
(MG-13) 

42 

106 

Lab  (AGREE) 

12. 

H-3113  Computers 

43 

2779 

Lab/Field 

13. 

Display  Console 

48 

314 

Lab/Field 

14. 

H-3182  Converters  D-A 

19 

102 

Lab/Field 

15. 

Magnetic  Tape  Units 

14 

50 

Lab/Field 

16. 

Infrared  Subsystem  in 

Fire  Control  Systems 
(MA-1) 

43 

99 

Lab  (AffiEF,) 

17. 

Rapid  Tune  Subsystem 

12 

41 

Lab  (AGREE) 

2.3.2  GOVERNMFWT  DATA  SOURCES 

The  search  for  and  acquisition  of  data  from  the  major  data  centers  main¬ 
tained  by  the  U.  S.  Government  Agencies  was  done  by  direct  contact.  Because 
of  the  vast  azoamt  of  data  collected  and  available  in  the  Government  data 
centers,  the  efforts  were  directed  to  collecting  only  adequate  and  useful 
data.  It  was  felt  that  a  visit  to  the  key  data  centers  was  a  necessity  to 
assure  a  thorough  analysis  of  the  type  of  data  collected  and  the  correct 
data  retrieval  methods. 

The  agencies  visited  were: 

1)  U.S.  Army  Maintenance  Command  Logistic  Data  Center  (USAMCLDC), 
Lexington,  Kentucky  (TAERS  data). 

2)  U.  S.  Navy,  Washington,  D.C.  ( 3M  data). 

3)  U.  S.  Navy,  Norfolk,  Virginia,  Statistical  Engineering  Branch, 
NAVSECNCRDIV. 

4)  U.  S.  Naval  Fleet  Missile  System  Analysis  and  Evaluation  Group 
(PMSAEG),  Corona,  California. 


5)  U.  S.  Air  Force  Logistics  Conmand,  Tinker  Air  Force  Base,  Oklahoma 
City,  Oklahoma  (66-1  data). 

There  was  no  reliability  data  available  os  electronic  equipment  at  the 
USAMC  Logistic  Data  Center  due  to  the  lack  of  control  of  serial  number 
assignment  to  end  items  and  usage.  Indies ting  devices  being  almost  non¬ 
existent  on  electronic  equipment  resulted  in  no  usage  time;  therefore,  no 
MTBF  data  is  available  threxigh  TAERS. 

An  additional  request,  suggested  by  the  USAMCLDC,  was  Bade  to  the  U.  S. 
Army  Electronics  Command  in  Fort  Monmouth,  Sew  Jersey,  Applications  Engineer¬ 
ing  Branch  which  resulted  in  MTBF  data f on  radio  «ets  but  was  unsatisfactory 
(n  too  small)  for  use  in  establishing  a  priori  distributions. 

The  Investigation  of  the  3M  data  collection  system  maintained  by  the 
Navy  yielded  no  reliability  data;  i.s.,  MTBF's.  The  Statistical  Engineering 
Branch  in  Norfolk.  Virginia  had  some  reliability  data  but  not  of  the  quality 
of  data  necessary  for  this  study.  A  trip  to  the  U.  S.  Naval  Fleet  Missile 
Systems,  Analysis  and  Evaluation  Group  in  Corona,  California  was  made  in 
search  of  MTBF  data  on  surface  missile  systems.  Again,  the  quality  (n  too 
small)  of  data  was  unsatisfactory  for  the  prior  distributions  requirement. 

A  visit  to  the  Reliability  Branch  at  Tinker  Air  Force  Base  in  Oklahoma 
City,  Oklahoma  resulted  in  acquisition  of  reliability  data  on  ground  communi¬ 
cations  radars.  Table  2, 3 *2.1  is  a  list  of  equipments  from  which  66-1  fora 
data  was  obtained.  All  of  this  data  was  based  on  a  f  1  d  time  T  of  4320 
hours. 


TABIE  2. 3. 2.1  TINKER  AIR  FORCE  BASE  ACQUIRED  DATA 


Assigned 
Data  Set 
Number 

Equipment 

n 

1. 

MTI  Reflector 

4l 

173 

2. 

HVPS  (High  Voltage 

Power  Supply) 

74 

475 

3- 

Synchronizer 

59 

265 

4. 

Oscilloscope 

51 

159 

5- 

Video  Amplifier 

4o 

92 

6. 

Synchronous  Power  Supply 

55 

183 

7. 

Search  Indicator 

58 

266 

a. 

Servo  Amplifier 

50 

194 
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2.3.3  DATA  SOURCES 


1.  Interceptor  Improvement  Progra 
Report,  May  196^,  Hughes  Airct 


IR  Reliability  Testing,  Final 
l^t  Company" 


2.  MMST  Production  Reliability  Testing,  "First  Reliability  Sampling 
Teat  foepcrt,"  July  1££>7,  “Second  Reliability  Sampling  Test  Report, 
October  1967,  "Fourth  and  Final  Reliability  Sampling  Test  Report," 
December  1967?  Hughes  Aircraft  Company. 

3.  MA-l/AR/ASQ-2?/M6-13  Interceptor  Improvement  Program,  "Rapid  Tune 
Reliability  Testing,1*"  November -De ceiSier  1965,  Hughes  Aircraft 
Company. 

V.  Iraa  Production  Reliability  Tests,  Final  Test  Report,  June  1968, 
Hughes  Aircraft  Company. 

2.3.^  SUMMARY  OF  DATA  AVAIIABILITY 

It  is  well-known  in  statistical  analysis  that  data  can  arise  in  two  ways 
i)  after  the  fact,  il)  as  a  result  of  an  experiment  designed  expressly  to 
answer  specific  questions;  for  example,  to  fit  a  prior  distribution.  The  re¬ 
sults  of  this  study  make  it  clear  what  to  do  if  data  is  to  be  gathered  ex¬ 
pressly  to  fit  prior  distributions.  However,  if  time  and  money  are  not 
available  far  the  second  approach  above,  already  existing  data  must  be  used. 
The  results  of  the  data  search  indicate  there  is  not  much  suitable  data 
available.  The  reasons  are  primarily  two 

1)  Not  much  data  is  available  on  large  numbers  of  identical  equipments. 

2)  Most  data,  already  in  existence,  involves  different  numbers  of 
failures  or  different  (fixed)  test  time  on  each  equipment. 

The  two  reasons  above  are  somewhat  different  in  character.  The  first, 
primer! ly,  causes  poor  fits  while  the  second  makes  it  difficult  to  apply  the 
methods  of  this  report  at  all. 

Since  large  amounts  of  data  cure  not  available  through  government  sources 
the  data  for  prior  fits  must  corns  from  Industry.  The  picture  may  appear 
unnecessarily  bleak:  the  results  of  Sections  b.O  and  5.0  indicate  data 
requirements  which  are  not  too  stringent. 
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SECTION  3.0  FUNCTIONAL  IZVEL  -  THE  CONDITIONAL  DISTRIBUTION 


In  using  Bayesian  methods  the  prior  distribution  receives  a  good  deal 
of  attention  -  and  rightly  so.  However,  in  any  analysis  involving  actual 
data  (X)  the  conditional  distribution  of  (X  !$)  must  also  be  known.  This 
conditional  p.d.f.,  sometimes  called  a  sampling  distribution,  receives  less 
attention  because,  usually,  more  is  known  about  it  than  about  the  prior 
distribution.  For  example,  if  several  machines  are  turning  out  a  large  number 
of  bolts  which  can  be  classified  only  as  good  or  bad  then  for  fixed  fraction 
defective  p  the  sampling  distribution  or  conditional  distribution  is  hyper- 
geometric  or  binomial  depending  on  the  finiteness  of  the  outputs.  Thus, 
often  the  physical  process  dictates  the  sampling  distribution.  In  other  cases, 
experience  has  dictated  what  to  expect  for  a  conditional  distribution.  That 
is  the  case  here.  We  have  assumed  that  the  important  Bayes  testing  applica¬ 
tions  will  be  made  at  the  equipment  (computer,  radar,  oscilloscope,  radio, 
etc.)  or  system  level.  In  this  case,  it  has  been  demonstrated  many  times 
both  by  limit  theorems  and  empirical  studies  that  if  the  random  variable  is 
time  to  failure,  a  good  descriptor  of  the  conditional  distribution  is  the 
exponential,  i.e. , 


f(x|*)=l/j  e  x,  8  > 0 

=  0  elsewhere. 

Then  the  distribution  of  the  number  of  failures  in  fixed  time  T  is 
Poisson,  i.e., 

-Ve  x 

f(x|ff)  =  e  (T/a)  t,  8>0 

x! 

x  m  0f  1,  ... 

-  0  elsewhere. 

It  should  be  understood  that  the  methods  given  here  are  illustrated  and  fully 
developed  for  the  above  two  conditional  distributions  but  that  in  general  the 
methods  are  applicable  to  any  conditional  distribution  provided  the  identifi- 
abillty  criterion  is  satisfied  (see  Appendix,  Section  9.4). 

If  it  should  turn  out  that  the  conditional  distribution  is  Weibuil,  then 
the  same  general  methods  would  apply  with  the  exception  that  the  prior 
distribution  would  be  two  (2)  variate.  The  major  change  then  would  be  that 
the  arithmetic  would  become  much  more  intricate. 
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SECTION  4.0  METHODS  OF  FITTING  ERIOR  DISTRIBUTIONS  AND  RESULTS 


SUMARY 

Objective 

The  Objective  of  this  section  is  to  develop  methods  of  fitting  prior 
distributions  to  8  *  MTBF  when  data  occurs  in  any  one  of  three  (3)  forms. 

Type  1  -  The  recorded  data  on  a  piece  of  equipment  is  number  of  failures 

occurring  in  a  fixed  time  T,  T  identical  for  all  s  equipments. 

Type  as  -  The  recorded  data  on  a  piece  of  equipment  is  the  observed 

MTBF,  6,  and  the  number  of  failures  on  each  Identical  equipment 

is  a  fixed  K. 

Type  2b  -  The  recorded  data  is  the  same  as  Type  2a  but  the  number  of 
failures  on  each  Identical  equipment  may  vary. 

Two  cases  were  considered:  the  family  of  the  prior  distribution  unspecified 
and  specified.  These  are  discussed  in  more  detail  in  Section  4.2.1 

Results 


Family  of  prior  distribution  unspecified,  all  types  of  data. 

Far  the  family  completely  unspecified  the  methods  found  were 
unsatisfactory  far  two  reasons,  the  first  being  the  most  important. 

1)  They  give  only  an  empirical  distribution  far  the  prior 
and  do  not  identify  a  family. 

2)  They  are  extremely  complex  to  apply. 

If  one  is  willing  to  assume  the  prior  distribution  belongs  to  the 
Pearson  family,  the  Krutchkoff -Rutherford  method  is  easy  to  apply 
but  requires  an  enormous  number  of  identical  equipment  for  good 
results. 

family  of  prior  distribution  specified,  all  types  of  data. 

When  the  family  la  specified,  the  steps  in  fitting  the  prior  dlstrl 
butlon  are  simple  enough. 

l)  Derive  the  marginal  distribution  undarths  assumed  prior  family. 


4.1  FITTING  THE  PRIOR  DISTRIBUTION  WHEN  THE  FAMILY  IS  UNSPECIFIED 

The  case  when  the  analyst  Is  unwilling  to  assume  any  knowledge  of  the 
prior  distribution  is  extremely  difficult  and  remains,  at  this  writing, 
essentially  an  unsolved  problem.  An  idea  why  can  be  obtained  by  looking  at 
the  classical  test.  Even  in  this  te3t  the  family  of  the  fitted  distribu¬ 
tion  must  be  specified.  The  methods  we  will  discuss  here  all  have  one  point 
in  commons 

/  .  -A  A  . 

They  all  use  the  data  (x^,  ...,  xq)  or  (0^,  ...,  0n)  to  estimate  the 

MOMENTS  of  the  prior  'distribution  and  then  apply  certain  uniqueness  theorems 
between  a  sequence  of  moments  and  a  distribution. 

The  extent  to  which  a  sequence  of  moments  determines  a  unique  probability 
distribution  AND  how  one  is  to  discover,  having  a  sequence  of  moments,  to 
which  probability  distribution  it  corresponds  is  a  difficult  problem  (see 
Shohat,  J.  A.  and  Tamarkin,  J.  D.  (1943),  The  Problem  of  Moments,  Amer.  Math. 
Society,  New  York).  The  methods  to  be  discussed  assume  some  of  the  difficul¬ 
ties  away. 

One  method  of  interest  is  given  by  John  E.  Rolph  (Bayesian  Estimation 
of  Mixing  Distributions,  Annals  of  Math.  Statistics,  Vol.  39>  No.  4,  August 
1968).  In  this  paper  Dr.  Rolph  assumes  a  prior  distribution  (uniform)  on  the 
family  of  prior  distributions.  Moreover,  he  assumes 

i)  The  parameter  space  0  is  limited  to  Lo,l],  i.e.,  0  i  0  k  1. 

ii)  The  conditional  distribution  f(x|0)  must  be  a  polynomial  in  0  and 
be  a  discrete  distribution. 

The  first  assumption  is  needed  to  use  the  theorem  that  a  probability 
distribution  on  [0,1]  is  uniquely  defined  by  its  moments.  The  first  assump¬ 
tion  also  permits  the  prior  distribution  on  the  family  of  prior  distributions 
to  be  assigned  to  the  moment  sequences.  The  second  assumption  permits  the 
marginal  distribution 

f(x)  =  y*f(x  |0)g(0)d0 

o 

to  be  written 

f(x)  »  E  axiui 
i 

wheie  is  the  i^^  moment  of  the  prior  distribution.  Thus,  a  sample 
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(xl>  *'*>  XQ)  froni  f(x)  can  bo  used  to  estimate  the  moments  of  the  prior  dis¬ 
tribution  and  hence  (because  of  uniqueness)  estimate  the  prior  distribution. 

The  estimate  constructed  in  this  manner  is  consistent.*  The  required  computa¬ 
tions  are  quite  Involved,  surely  requiring  a  computer  program  since  a  number 
of  relatively  high  order  determinants  are  involved.  This  method  is  not 
directly  useful  here  since  the  restriction  of  0  to  Co,l]  is  untenable.  How¬ 
ever,  it  is  mentioned  since,  in  the  exponential  case,  the  reliability  function 

R(T)  =  e~T/0  T  fixed, 

always  satisfies  0  ^  R(T)  =*  1.  Thus,  a  prior  distribution  could  be  fitted 
to  R(t)  and  then  a  change  of  variable 

9  -  -T/inR 

leads  to  a  prior  distribution  on  6-  It  is  not  possible  to  do  this  in  this 
study  because  repeated  observations  with  T  fixed  are  not  available  from  the 
data  search.  It  should  be  noted  though  that  if  it  were  possible  to  make 
successive  Bernoulli  trials  then  the  conditional  distribution  of  f(x|R)  Is 
a  polynomial  in  R  and  the  method  (with  the  change  of  variable  Q  =  -T/inR) 
could  be  used  to  fit  g(0). 

Howard  G.  Tucker  has  given  a  method  for  estimating  the  prior  distribution 
when  the  conditional  distribution  f(x|6)  is  Poisson  (An  Estimate  of  the  Com¬ 
pounding  Distribution  of  a  Compound  Poisson  Distribution,  Theor.  Prob.  Appl. 

8,  195-200,  1963).  Here,  the  observed  random  variable  x  is  not  lifetime  but 
the  number  of  failures  occurring  in  fixed  time  T.  Thus,  it  is  assumed  that 
random  samples  (x]_,  ...,  Xq)  are  available  from  the  marginal  (discrete)  dis¬ 
tribution  f(x).  Again,  the  uniqueness  of  the  moment  sequence  for  the  prior 
distribution  is  used  to  obtain  a  consistent  estimate  of  the  prior  distribution. 


Nov,  we  assume  throughout  this  study  that  for  lifetimes,  the  conditional 
distribution  of  lifetime  given  0  is  known  and  that  though  each  of  the  n  0's 
is  unobservable  (because  6  for  each  equipment  is  unknown)  samples  are  available 
from  f(x|0),  where  X  represents  lifetime.  The  sample  size  from  f(x|0)  will 
be  denoted  by  for  the  i^h  0.  We  sure  particularly  concerned  with  methods 


which  can  deal  with  the  conditional  distribution 


( V  •*) 


since  often  the 


only  information  available  to  the  analyst  ie  the  pair  (§,K. )  rather  them  the 
lifetimes  (x^ ,  ...,  x^). 


The  case  large  i»l,  ...,  n(e.g.,  >  30  all  i)  can  be  dismissed  more 

or  less  out  of  hand  since  (see  Section  9-5)  the  Joint  (marginal) 
distribution  f (^i,  • • • y  ^n)  can  be  used  directly  to  fit  the  prior 

distribution  g(0).  Before  proceeding  to  other  cases,  it  is  well  to 


Roughly,  a  consistent  estimate  means  that  in  the  limit,  i.e.,  as  n  -•00 
the  true  prior  distribution  will  be  known  exactly. 
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c^.ear  up  some  terminology.  In  reliability  work  the  estimation  of  the  prior 
distribution  is  called  fitting  a  prior  distribution.  In  some  of  the  more 
theoretical  work,  it  is  called:  estimating  the  mixing  distribution  or  some¬ 
times:  estimating  the  compounding  distribution.  Thus,  in  the  statistical 
literature,  the  prior  distribution  is  sometimes  called  mixing  distribution  or 
the  compounding  distribution. 

The  final  method  we  discuss  for  fitting  an  unspecified  prior  distribu¬ 
tion  is  the  Krutchkoff -Rutherford  (hereafter,  K-R)  method.  Strictly  speaking, 
the  K-R  method  is  not  usable  when  the  prior  distribution  is  completely  un¬ 
specified.  The  K-R  method  requires  that  the  prior  distribution  belong  to 
the  Pearson  family  of  curves.  This  is  not  very  restrictive  since  the  Pearson 
family  is  a  rather  large  family  (e.g.,  it  includes  the  inverted  gamma  distri¬ 
bution).  It  turns  out  that  members  of  the  Pearson  family  are  uniquely  deter¬ 
mined  by  their  first  four  moments.  Thus,  computations  are  greatly  simplified 
(still  requiring  a  computer  program  though)  as  regards  the  previous  two 
methods  discussed.  As  in  the  other  methods,  the  K-R  method  leads  to  consis¬ 
tent  estimates.  There  are  two  serious  shortcomings  with  the  K-R  method: 

i)  It  is  not  always  usable. 

ii)  It  has  large  and  unknown  sampling  errors. 

The  second  shortcoming  is  common  to  all  the  methods  discussed  but  is 
particularly  true  of  the  K-R  method  because  it  uses  only  the  first  four 
moments.  The  first  shortcoming  above  is  the  most  serious.  If  the  sample 
second  central  moment  is  negative  the  K-R  method  cannot  be  applied.  The  K-R 
method  was  tried  on  seventeen  sets  (seventeen  different  types  of  equipment) 
of  data.  Only  four  sets  resulted  in  an  estimate  of  the  prior  distribution, 
the  other  thirteen  resulting  in  negative  sample  second  moments.  The  K-R 
method  is  discussed  in  detail  in  Section  4.3*1  of  this  report. 

4.2  METHODS  OF  FITTING  PRIOR  DISTRIBUTIONS  WHEN  THE  FAMILY  IS  SPECIFIED 
4.2.1  INTRODUCTION 

The  problem  of  fitting  a  prior  distribution  is  much  simpler  when  the 
prior  distribution  family  is  specified  than  when  it  is  unspecified.  Once  a 
distribution  family  (e.g.,  inverted  gamma),  is  assumed  for  the  prior  then 
fitting  the  prior  is  a  two-fold  procedure: 

1)  Obtaining  estimates  of  the  parameters  of  the  prior  distribution. 

2)  Testing  the  validity  of  the  assumed  prior  by  a  goodness -of -fit 
test. 

Before  discussing  how  l)  and  2)  are  carried  out  the  general  setup  is 
recalled: 
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There  is  a  prior  distribution  g(0)  on  the  parameter  0,  where  0  is 
always  considered  the  mean  of  an  exponential  time -to -failure  distribution 

(i.e.,  f(x|0)  **  (1/0)  e  ®).  Sample  value s  from  some  observable  random 
variable  are  available.  In  this  study,  only  two  types  of  observable 
random  variables  are  of  sufficient  practical  interest  to  be  investigated: 


Type  1.  The  observed  random  variable  is  the  number,  X,  of  failures  of 
a  unit  occurring  in  a  fixed  time  T.  X  is  a  discrete  variable, 
taking  on  only  the  values  0,1,2,  ...  .  Observations  on  X  are 
obtained  by  pitting  n  units  on  test  for  time  T,  and  recording 
the  number  of  failures  for  each  unit  (x^,  ...,  x^} . 

Type  2.  The  observed  random  variable  is  the  sample  MTBF,  0.  Two  cases 
are  possible  here: 


Type  2a.  Each  value  of  0  i^  based  on  a  fixed  number  of  failures,  K. 
Sample  values  of  0  are  obtained  by  putting  n  units  on 
test  and  recording  K  failure  times  x.  . . . ,  x .  v  for 

JL  j  -L  1  y 

each  unit  i  =  l,  ...,  n,  and  then  computing 


(Note  that  since  y^is  a 


,  X.  v) ,  one  has  no 
%  A 


sufficient  statistic  for  (X.  , , 

A  1>1  . 

need  of  the  x.  *s  once  0.  has  been  computed).  0  is  a 
y  <J  1 

continuous  random  variable  defined  on  (assuming, 

of  course,  that  0  also  is). 


Type  2b. 


The  values  of  0  are  based  on  different  values  of  K.  This 
is  the  case  of  practical  interest.  Data  in  the  form  of 
Case  2a  does  not  occur  in  practice,  but  Case  2a  is  in¬ 
cluded  be  anise  its  treatment  is  analytically  simp!:-:  and 
sheds  light  on  Case  2b.  Foj^  this  case,  the  observed  data 

is  in  the  form  of  n  pairs  (0. ,K. ),  where  K.  is  the  number 

1  1  1  A 

of  failure  tines  used  for  the  computation  of  0^. 


Whether  the  observed  data  is  of  Type  1  or  Type  2, 
the  situation  is  the  same:  the  observed  data  is  not  from 
the  prior  g(0),  tut  from  the  marginal  distribution  f(x) 

(or  f (0 ) ) .  Hence,  the  data  cannot  be  fitted  "directly’’ 
bo  g(e) ,  However,  when  g(0)  is  specified,  f(x)(f(e))  is 
uniquely  specified,  since 


00 

f(*)  «  j g(0)f(«|0)d0 
0 
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and  f(x|e)(f(£|e))  is  specified  as  a  particular  Poisson 
(gamma)  distribution  by  the  exponent  Dial  assumption. 

Since  the  parameters  of  f(x)(f(o))  include  the  parameters 
of  g(9 )f  one  need  only  estimate  the  parameters  of  the 
marginal  to  obtain  estimates  of  the  prior  parameters. 
Also,  because  of  the  unique  association  between  the 
prior  and  the  marginal,  a  goodness  of  fit  test  of  the 
marginal  also  serves  as  a  goodness  of  fit  test  for  the 
prior-. 

In  summary,  a  specified  prior  yields  a  specified  marginal  contain¬ 
ing  the  same  parameters.  Since  the  observed  data  comes  from  the 
marginal,  one  can  fit  the  data  to  the  marginal  using  any  appropriate 
classical  methods  of  parameter  estimation  and  goodness  of  fit  tests. 

In  the  process,  the  prior  g(e)  is  also  fitted. 


In  the  following  sections,  specific  formulas  are  derived  for  fitting 
the  prior  when  the  specified  prior  family  is  the  inverted  gamma  distri¬ 
bution 


g(e)* 


X.4  i 


A*h) 


(4. 2. 1.1) 


The  inverted  gamma  prior  distribution  is  selected  because  it  is 
a  two  parameter  distribution,  very  flexible,  and  the  natural  conjugate 
when  the  conditional  distribution  is  exponential.  Also,  when  the 
observed  data  is  of  Type  1  or  Type  2  -above,  the  marginal  density  is 
available  in  closed  form. 

The  method  selected  for  estimating  parameters  is  the  method  jf 
matching  moments.  For  goodness  of  fit,  the  test  is  used. 


4.2.2  PARAMETER  ESTIMATION 

4. 2. 2.1  PARAMETER  ESTIMATION  FOR  TYPE  1  DATA 

When  the  observed  statistic  is  the  number  of  failures,  X,  occurring 
in  fixed  time  T,  then  its  distribution  is  Poisson, 

f(*le)  .  e~T/9(T/l£  T,e  >  o  (4.2.2.1.1) 

X  • 

x  c  0,1,  • •  • 

since  the  conditional  distribution  of  time  to  failure  for  fixed  6  is 
exponential.  With  a  specified  inverted  gamma  prior  distribution 
■  or  ;ulr  4. 2. 1.1),  the  marginal  distribution  of  X,  derived  in  Section 

;  ■  is  given  by 
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f(x) 


p 


x  -  0,1,  ...  (4. 2. 2. 1.2) 


i 


I 


This  is  a  negative  binomial  distribution  with 


Mean  >•  E<X)  -  2? 

(4.2.2.1.3) 

Variance  -  -  “gisl 

(4. 2. 2. 1.4) 

ECX2)  -  (E(X))2  ♦  nx2  -  Mjpg)  ♦  ^ 

(4.2.2.1.5) 

Suppose  n  sample  values  of  X,  (r.,  ...,  x^)  are  observed. 

Then  b;T  the  method  of  matching  moments,  we  equate 

(Vxjj/n  -  tf/« 

(Eh>- (ff)* 

(4.2.2.1.6) 

(4. 2. 2.1.7) 

The  values  of  «  and  X  which  satisfy  Formulas  4. 2. 2. 1^6 
4. 2. 2. 1-7  are  the  required  estimates,  and  are  denoted  by  af 
The  solution  is 

A  —  T 

a  •  "  5" '  — :g . 

IS)  . 

n  \  n  /  n 

r 

(4. 2. 2. 1.8) 

(4.2.2.1.9) 

These  formulas  trill  be  used  in  Section  4.3  to  obtain  esttmte* 
of  3  and  A.  far  eight  sets  of  field  data. 


4. 2.2.2  PARAMETER  tirgmkTXXM  F3B  TTPE  2s  l ATk 

A  Type  2s  data  occurs  v&en  the  observed  statistic  is  the  Mugpla  WB?, 

8,  where  the  observed  number  of  failure*  is  a  fixed  integer  K.  The  cooiltiond 

distribution  it  gsnaa,  i.e., 


*K(ele) 


k\K  1  AK4  T 

e)  fW  (e) 


(4. 2. 2.2.1) 


because  of  the  assumption  of  an  exponential  distribution  of  time  to  failures. 

With  a  specified  inverted  gsuan*.  prior  distribution  (Formula  4.2.11),  the 

A 

marginal  distribution  of  0,  derived  in  Section  9 -.2,  is  given  by 

(^)-  Ko’  (*-2-£-2-2> 
This  is  an  inverted  Beta  distribution  with 


t  \  >  1 


(4. 2. 2. 2. 3) 


Kft  -  ¥ 


X  >  2 


k2.2.2.4) 


A  .A  A  x 

Suppose  n  sample  values  of  0,  (61>  • • • >  0Q),  are  observed,  each 

based  on  K  failures.  Then  the  method  of  matching  moments  leads  tc  the  pair 
of  equations 


i  a 
"n  “ 


(4. 2. 2. 2. 5) 


E'ef 


K+l  _ <•? 

k  TTOT 


(4. 2.2.2. 6) 


A  ^  A  K 

For  convenience,  we  eet  m.  »  i  ,  nu  ■  *  _ i 

n  «  R*1!  n 


•a  n, 

Then  the  solution  (ar,X)  is  given  by 


A  2 

X  “  1 C-TTJ 
“2"“! 


(4. 2. 2. 2.7) 
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a 


A  A  % 

a  -  rn^(x-l) 


(4. 2. 2, 2. 8) 


4. 2. 2. 3  PARAMETER  ESTIMATION  PCS  TYPE  2b  DATA 

A 

A  When  the  observed  statistic  is  d,  bat  the  number  of  failures  frees  which 
ft  is  computed  varies  from  unit  to  unit,  a  ■edification  of  the  Method  In 
Section  4. 2. 2. 2  can  be  used  to  obtain  estimates  of  a  and  A. 

Type  2b  data  is  available  in  the  fora  of  n  observations,  each  observation 
being  a  pair  (*0i,K1),  i  -  1,  n,  where  the  positive  integer  K ^  is  the 

number  of  failures  used  to  compute  e  .  Now,  far  any  value  of  K. ,  if  we 

A  1  1 

consider  the  random  variable  6  based  on  failures,  we  obtain  from  equations 

4. 2. 2.2. 3  end  4. 2- 2. 2. 4. 


,A. 

y®) 


{4.2-2. 3.1) 


(h  A 

\ V1  / 


(4. 2.2.3. 2) 


Since  the  above  formulas  are  valid  far  any  value  of  K^,  the 
following  procedure  for  estimating  a  and  \  from  the  sample  data 

^0l'K,)5i-i  seems  reasonable:  Solve  the  equations 


(4. 2. 2. 3*3) 


£  1  ft  2 


(4. 2.2. 3*4) 


for  cr  and  X-  For  convenience,  we  set  m^ 


/ 
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A 

“2 


X  / 

£  /  D*  T5sen  solution  (a,£)  is  given  by 


A  A  ‘ 
A  ^"“l 

X  *  7T7TT 

■S'*3! 

a  *  m,  (£-1) 


(4. 2. 2. 3*5) 

(4. 2=2. 3*6) 


Bote,  that  ve  have  obtained  point  estimates  of  a  and  X  without  haring 
to  consider  what  the  form  cf  the  marginal  density  f(§)  is  when  the  Ki's 
are  different.  This  was  possible  because  the  expectations  in  Formulas 
4. 2. 2. 3*1  and  4. 2. 2. 3*2  are  independent  of  K^*  However,  when  goodness  of 
fit  is  discussed  in  the  next  section,  it  will  become  necessary  to  specify 
a  model  for  the  marginal  distribution. 

4.2.3  GOODNESS  OF  FIT 

4.2.3-1  THE  X2  TEST 
2 

The  X  test  was  used  in  this  study  to  test  the  goodness  of  fit  of  the 
estimated  prior  distributions.  This  section,  will  outline  the  steps  involved 

in  taking  the  X?  test.  The  theory  behind  the  X2  test  (anL  also  the  method 
of  matching  momenta  used  in  Section  4.2.2)  can  be  found  in  many  standard 
textbooks  on  Statistics. 

Suppose  we  have  specified  a  family  of  probability  distributions  (e.g., 
negative  binomial  of  form  (4. 2- 2. 1.2)),  and  have  used  sample  data  to  estimate 
q  unknown  parameters  of  the  distribution.  To  test  goodness  of  fit: 

1)  Select  a  significance  level  p  (e.g.,  p  »  .90,  P  *  *95)* 

2)  Divide  the  range  of  the  data  into  cells,  so  that  at  least  5  sample 
values  lie  in  each  of  the  c  cells. 

3)  Count  the  number  of  observations  in  each  cell,  and  compute  the 
expected  number  of  observations  in  each  cell  under  the  hypothesis 
that  the  data  arose  from  the  distribution  with  the  estimated 
parameters. 

4)  Denoting  the  observed  and  expected  frequencies  In  each  cell  by  0., 
E^,  1*1,  ...,  c,  compute 
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(4. 2.3. 1.1) 


X 


c  (o . -e. y 

E  -V- 

Ei 


2  f  w  p 

5)  Compare  k  to  the  p  quantile  of  the  y  distribution 

2 

with  c-q-1  degrees  of  freedom  (denoted  X  p,c-q-l). 


If 


v2  2 
x  >  X 


p,c-q-l 


(4. 2.3. 1-2) 


reject  the  hypothesized  distribution  family  (i.e.,  a  bad  fit 
has  been  demonstrated). 


If 


1(2  <  (4. 2.3. 1*3) 

then  accept  (i.e.,  good  fit). 

We  now  return  to  the  case  of  fitting  the  inverted  garaa  prior  for 
different  types  of  marginal  distributions.  In  each  case,  the  X2  test  is 
taken  against  the  marginal,  since  that  is  the  distribution  from  which  the 
sample  data  arises. 

4. 2. 3-  2  THE  X2  TEST  FOR  TYPE  1  DATA 

When  the  observed  statistic  is  the  number  of  failures,  X,  occurring  in 
fixed  time  T,  Formula  (4. 2. 2. 1.2)  is  used  to  compute  f(i)  -  P(X*i)  for 
i  ■  0,1,2,  ...  .  For  a  sample  of  size  n,  the  expectation  of  X»i  is  nf(i). 

Once  the  integers  in  the  range  of  X  are  divided  up  into  cells,  the  expected 
number  of  observations  in  each  cell  is  easily  summed. 

With  c  cells,  the  x  test  taken  with  (c-3)  degrees  of  freedom, 
since  2  unknown  parameters  (a  and  X)  were  estimated. 

4. 2. 3-  3  THE  X2  TEST  FOR  TYJE  2a  DATA 

When  the  observed  statistic  is  e  based  on  a  fixed  number,  K,  of  failures, 
then  the  marginal  distribution  fK(@)is  the  inverted  Beta  distribution  given 

in  Formula  4. 2*2. 2. 2.  For  a  sample  size  of  n,  the  expected  number  of  obser¬ 
vations  in  a  x2  cell  with  upper  and  lower  end  points,  U  and  L,  respectively, 
18 

n  J  fK(e)de  (4. 2.3. 3.1) 

L 
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The  integral  is  not  available  in  closed  form,  but  can  be  evaluated  by 
raking  a  transformation  to  the  Beta  distribution  and  using  a  table  of  incom¬ 
plete  Beta  functions.  If  the  tables  do  not  have  the  appropriate  parameter 
values  for  a  particular  case,  it  is  always  possible  to  use  the  computer  tech¬ 
niques  of  either  numerical  integration  or  simulation  of  the  distribution. 

2 

Since  only  2  unknown  parameters  are  estimated  (K  is  known),  the  X  test 
is  taken  with  degrees  of  freedom  3  less  than  the  number  of  cells. 

4.2. 3. 4  THE  X2  TEST  FOR  TYPE  2b  DATA 

A  A  A 

When  the  observed  data  consists  of  n  pairs  jXfl. ,K . )},  and  a  and  X  are 

*  *  /\ 

estimated  as  in  Section  4. 2. 2. 3,  a  modgl  for  the  marginal  distribution  f(y) 
must  be  specified  in  order  to  take  a  X  te^t.  This  model  must  account  for 
the  different  K^'s  on  which  the  values  of  0  are  based.  The  most  reasonable 

model  seems  to  be  the  mixed  population  model,  in  which  the  marginal  density 
takes  the  form 


ftit1)(e)  •  piV8>  *  -  *  ptye)’ 


(4.2.3.11.1) 


where  t  is  the  number  of  distinct  K. 'e,  and  p,  is  the  probability  that  a 

A  1  1  t 

random  observation  of  0  is  based  on  K.  failures  (note  £  Pi  *  l) •  The 

i*l 

f  (©) *  s  are,  of  course,  the  inverted  Beta  marginal  densities  (Formula 
1  A 

4. 2. 2. 2. 2)  for  the  case  of  6  being  based  on  failures.  The  mixed  population 
model  assigns  relative  "weights"  to  the  "single  population"  marginals  f  (e), 

1 

„  „  A 

in  accordance  with  the  prior  probabilities  p.  of  a  sample  6  being  based  on 
K.  failures.  It  is  easily  verified  that  t!ir  -i(0)  is  indeed  a  density: 

/  ■ ■  f  -  £%  /  y6>^ 


(4. 2. 3- 4. 2) 


The  model  contains  (t+2)  unknown  parameters  which  must  be  estimated 
2 

before  a  X  test  can  be  taken:  or,X,Pi,  •••,  p+.  The  estimation  of  a  and  X 

x  va 

has  already  been  discussed.  For  a  sample  of  n  0's,  the  pi  are  estimated  by 


(4.2.3  ^-3) 


A 

P« 


number  of  e's  computed  from  failures 


for  i*l ,  * • . ,  t  > 

2 

To  compute  the  expected  number  of  observations  In  a  X  cell  with  upper 
and  lower  end  points  U  and  L„  respectively,  we  note  that 


U  U 

fK  (e)de  ,  (k.2.3.k.k) 

L  1  1  L  1 


and  the  evaluation  of  the  integrals  on  the  right  was  discussed  in  Section 

4. 2. 3- 3- 


With 


c  \  cells,  the  x  t^st  is  taken  with  c-t-3  degrees  of  freedom. 


since  t+2  parameters  of  f „  ,(0)  sire  estimated. 

*-Ki- 


4.3  RESULTS  AND  DEVELOPMENT  OK  PRIOR  DISTRIBUTIONS  FOR  DATA  COLI£CTED 

4.3-1  RESULTS  OF  FITTING  THE  PRICK  DISTRIBUTION  WHEN  THE  FAMILX  IS 
UNSPECIFIED 


The  methods  for  fitting  the  prior  distribution  when  the  funily  is  un¬ 
specified  were  discussed  in  Section  4.1.  Only  one  of  these  methods  was  used 
for  fitting  actual  field  data  in  this  study:  the  Krutchkoff -Rutherford 
method.  The  necessary  type  of  data  for  using  the  K-R  method  (pairs 

was  available  for  all  17  sets  of  data  collects!  (See  Section  ?) .  Table 
4.3.1  summarize e  the  results  of  applying  the  K-R  method  to  the  17  sets  of 
data  (numbered  as  in  Section  2).  The  formulae  at  the  right  are  those  used  in 

calculations.  The  m. 'a  are  the  unbiased  estimates  of  the  moments  of  g(9) 

1  A 

that  the  K-R  method  calls  for.  The  are  estimates  of  the  central  moments 


of  g(9)  obtained  by  using  the  relationships  between  the  central  moments  and 
the  moments  about  0  that  held  for  any  distribution.  Finally,  the  estimates 


oi  the  u.  are  used  to  obtain  estimates  of  A 


,  #P2  and 


which  are  used  in 


accordance  witn  criteria  in  Eldsrton't*  book  (Reference  2)  to  Identify  uniquely 
the  appropriate  Pearson  curve  type. 


In  Table  4,3.1,  Pearson  curve  types  are  identified  for  only  Data  Set* 

Nos.  2,  U,  12,  and  14,  which  are  fitted  as  Type*  VI,  I,  I,  and  IV  respective¬ 
ly.  For  the  13  other  cases,  the  K-R  method  falls  to  work,  since  the  esti¬ 
mates  of  we  negative.  The  failure  to  fit  a  prior  in  13  out  of  17  cases, 

is,  of  coarse,  no  fault  of  the  data,  but  is  inherent  in  the  K-R  method.  Even 
in  the  four  cases  in  which  a  prior  distribution  is  identified,  the  results 
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are  not  very  credible  because  of  the  large  and  unknown  sampling  errors  in¬ 
herent  In  the  method.  To  see  Just  how  "bad"  the  errors  are,  some  simulations 
of  the  K-R  method  (to  be  discussed  in  Section  5)  were  carried  out.  The  re¬ 
sults  indicate  that  the  K-R  method  gives  hopelessly  bed  results  when  used  on 
data  of  any  practical  sample  size. 

A  list  of  the  computer  program  used  far  the  K-R  method  is  given  in 
Section  10. 

4.3.2  RESULTS  OF  FITTING  THE  HUCR  DISTRIBUTION  WHEN  THE  FAMILY  IS  SPECIFIED 

The  results  of  fitting  an  Inverted  gamma  prior  distribution  are  shown 
in  Table  4. 3. 2.1.  The  8  sets  of  field  data  used  are  those  given  in  Table 
2. 3. 2.1.  In  all  8  sets  of  data,  all  units  were  put  on  test  for  time  T  ■  4320 

hours  and  the  number  of  failures,  X,  recorded.  (None  of  the  other  9  sets  of 

data  were  based  on  fixed-time  testing).  Therefore,  the  method  of  parameter 
estimation  used  is  that  given  in  Section  4. 2. 2.1  (Type  1  data).  Table 

4. 3. 2.1  gives,  for  each  set  of  data,  the  sample  mean  £x  /n,  the  sample  variance 

p  ^  ^  * 

/n  -  fco^/n)4  ',  and  the  Inverted  gamma  parameter  estimates  a  and  \,  calcu¬ 
lated  by  Formulas  4. 2. 2. 1.5  ajjd  4. 2. 2. 1-9,  respectively.  The  last  4  columns 
in  the  table  pertain  to  the  X  goodness-o£-fit  tests  as  described  In  Sections 

4. 2. 3.1  and  4. 2. 3. 2.  TJje  number,  c,  of  X  cells  selected  is^given,  along  with 
the  computed  value  of  X  ,  and  a  desig.iation  of  whether  the  \  test  passes  or 
fails  (with  c-3  degrees  of  freedom)  at  both  significance  levels  p  *  .99  and 

P  *  .90. 

2  From  Table  4. 3. 2.1,  one  can  see  that  at  the  .99  significance  level,  the 
X  '  test  is  passed  in  7  out  of  6  cases.  Since  each  case  is  far  a  different 
type  of  equipment,  these  results  indicate  a  general  applicability  of  the 
inverted  gamma  prior  distribution  cn  0  for  a  large  rang*  of  equipment  types. 
Figures  4. 3. 2.1  through  4.3.2.14  give  the  empirical  and  theoretical  marginal 
distributions  and  the  theoretical  inverted  gamma  prior  distributions  which 
have  been  fitted. 

In  Section  10,  there  is  a  list  of  the  computer  program  used  to  calculate 
A  A  2 

a,  X,  and  X  from  the  field  data.  Included  is  a  list  of  the  raw  data  for 
each  of  the  8  cases,  indicating  the  grouping  into  X2  cells.  A  listing  of  the 
out nut  of  the  program  is  also  included. 


An  attempt  was  made  to  fit  a  set  of  field  data  to  a  Weibull  prior  dis¬ 
tribution: 


/  »  0  8-1  -0^® 
g(b)-  ~  6  e  '  • 


(4. 3. 2.1) 


flats  Set  No.  4,  with  estimated  inverted  gums  prior  parameters  a  ■  3609.6  and 

X  «  2.6 05,  was  selected.  Weibull  parameters  were  selected  to  yield  the 
sane  mean  and  variance  as  the  above  inverted  ga am  distribution.  The  esti¬ 
mates  of  the  Weibull  scale  and  shape  parameters  are  then,  respectively. 


3.2 


( 


a  *  385.65,  £  »  .7656.  In  Table  4. 3, 2. 2,  the  Xc  cells  for  Data  Set  Ho.  4  are 
shown,  giving  for  each  cell  the  observed  end  the  expected  number  of  observa¬ 
tions  under  both  the  inverted  gamma  hypothesis  and  the  Weibull  hypothesis. 

The  expected  values  for  the  inverted  gpunaa  hypothesis  are  computed  as 
described  in  Section  4.2. 3.2.  For  the  Weibull  case,  however,  the  marginal 
distribution  is  not  available  in  closed  fora,  so  it  had  to  be  simulated. 

For  T  *  4320,  a  -  385, 65,  and  0  *  .7856,  10,000  random  values  of  X  (number  of 
failures  in  time  T)  were  drawn  in  a  two-stage  process:  first,  a  random  0  is 
drawn  from  the  Weibull  prior,  then  a  random  number  of  failure  tines  in  T  given 
0  is  drawn.  The  relative  frequencies  in  10,000  trials  of  the  events 
1X-0 ,  X*l,  ...}  are  used  to  approximate  the  marginal  distribution  f(x). 

2 

The  calculated  values  of  X  are  1.315  end  16.482  far  the  inverted  gamma 
hypothesis  and  the  Weibull  hypothesis,  respectively.  Since  in  each  case  two 
unknown  parameters  are  estimated,  both  x  tests  are  taken  with  4-1-2  ■  1 
degree  of  freedom.  For  the  inverted  guana  case,  the  X2  test  is  passed  at 
significance  level  p  ■  .90,  and  even  with  p  as  low  as  it  is  passed  after  .70. 
The  Weibull  hypothesis,  however,  fails  the  X2  test  even  at  level  p  »  .9999. 
(Note  that  p  is  the  probability  of  acceptance  when  the  hypothesic  is  true, 
and  thus  it  is  easier  to  pass  a  X2  test  when  p  is  high).  Sinue  the  hypothe¬ 
sized  Weibull  prior  fails  the  x2  teat  so  decisively,  this  result  significantly 
helps  to  validate  the  assumption  of  an  inverted  prior. 


In  this  study  no  type  2a  data  was  uncovered  although  by  a  designed  teat 
such  data  certainly  could  be  gathered.  Much  of  the  data  gathered  in  this 
study  was  of  type  2b ;  in  fact,  seventeen  (17)  sets.  Thus,  the  mixed  marginal 
distribution  model  had  to  be  applied.  This  model  was  discussed  in  Subsection 
4. 2. 3. 4.  In  order  to  apply  the  X2  test  for  type  2b  it  is  necessary  to  have 
the  number  of  cells  in  the  X2  test,  c,  such  that  c  *  t+4.  In  none  of  the 
seventeen  (17)  data  sets  was  the  sample  size  n  large  enough  so  that  at  least 
t+4  cells  were  available.  This  is  not  as  severe  a  limitatio.'  as  it  first 
might  appear.  The  number  of  identical  equipments  n  should  be  as  large  as 
possible  for  good  fits  and  if  n  is  relatively  large  with  respect  to  the  number 
of  distinct  K's,  i.e.,  t  the  mixed  model  can  be  applied. 


4.3.3  REMARKS  ON  TEE  FITTED  PRICK  DISTRIBUTIONS 


ifcta  sets  1,  2,  4,  5,  6,  7,  3  have  been  shown  to  be  well  described  by 
an  inverted  gamma  prior  distribution.  These  seven  data  sets  le present  seven 
types  of  equijment.  The  question  arises  as  to  the  decree  of  appli¬ 
cability  of  these  results  to  other  similar  equipments.  For  example,  Data  eet 
4  represented  an  oscillosc There  are  many  different  types  of  oscilloscopes 
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PaiOB  DISTRIBUTIONS 


Inverted  for  Weibull 


by  various  manufacturers  and  one  may  ask:  v.ould  taese  oscilloscopes  have 
the  same  prior  distribution  as  that  of  Data  set  4?  (and  hence,  the  same 
estimated  prior  distribution  as  we  have  achieved  could  b?  used).  The  answer 
is  strictly  speaking,  no.  However,  his  is  probably  speaking  too  strictly 
for  practical  purposes.  It  is  very  likely  that  equipments  similar  to  those 
used  in  this  report  will  have  piuor  distributions  of  the  same  family  (i.e., 
inverted  gamma)  and  hence,  one  need  only  estimate  the  parameters  and  not  go 
through  the  sometimes  tedious  process  of  fitting  the  prior  distribution  itself 
It  may  also  turn  out  that  even  if  the  prior  distributions  do  differ,  they  may 
be  combined  into  one.  This  is  actually  not  recommended  because  the  combina¬ 
tion  of  two  inverted  gamma  distributions  in  this  way  leads  to  a  mixed  distri¬ 
bution  which  it  not  inverted  gamma.  C '.early,  there  are  two  important  areas 
for  future  study 

i)  The  applicability  of  fitted  prior  distributions  to 
similar  equipments. 

ii)  The  feasibility  of  combining  equipments  into  one 
prior  distribution. 
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SECTICN  5.0  REQUIREMENTS  FOR  SUITABIE  A  PRIORI  DATA 


.  UMMARY 
Objective 

The  objective  of  this  section  is  to  develop  requirements,  qualitative 
and  quantitative,  on  the  type  and  amount  of  data  suitable  for  fitting  a 
prior  distribution  to  e  »  MTBF. 

Results 


Regarding  type  of  data,  it  was  found  that  virtually  all  types  of  failure 
data  are  suitable  if  there  is  enough  of  it.  Even  observations  on  probability 
of  survival  can  be  used.  The  "preferred"  types  are  the  aforementioned  type  1 
and  type  2a  but  type  2b  can  also  be  used.  The  results  for  the  amount  of  data 
are  summarized  below. 

Type  1  -  number  of  failures  in  fixed  time  T. 

n  *  30,  T  large  enough  to  obtain  variation  in  the  number  of 
failures  observed  on  each  equipment.  For  example,  T  large  enough 
so  that  the  number  of  failures  are  not  all  the  same  for  all  n  equipments. 

Type  2a  -  observed  MTBF,  ti,  K  identical  for  all  n  equipments. 

n  >  20,  K  >  10. 

A 

Type  2b  -  observed  MTBF,  a,  K  not  the  same  for  all  equipments. 

The  number  of  identical  equipments  n  should  be  large  enough  so 
that  n  *  5t  +  20  where  t  is  the  number  of  distinct  K^'s. 


48 


Along  with  the  development  of  methods  of  fitting  prior  distributions, 
an  important  goal  of  this  study  is  to  establish  minimal  data  requirements 
for  using  these  methods.  There  are  two  aspects  of  the  suitability  of  data  to 
be  considered: 

i)  The  type  of  data. 

ii)  The  amount  of  data. 

5.1  THE  TYPE  OF  DATA 

Failure  data  can  occur  in  several  forms .  We  divide  these  into  two 
classes:  attributes  and  variables  data.  The  attributes  situation  occurs 
when  an  equipment  is  operated  and  its  survival  or  nonsurvival  for  time  T 
(usually  mission  time)  is  observed.  Then  usually  the  binomial  distribution 
describes  such  observations  and  the  parameter  probability  of  survival  in  the 
binomial  distribution  is  the  reliability  function 

R(T)  -  e'T/0,  T  fixed.  (5*1.1) 

Now,  since  the  MTBF  »  6  is  being  considered  a  random  variable  in  this  study, 
then  so  is  R(t)  and  either  Rolph's  method  or  the  K-R  method  (See  Section  4.1) 
can  be  used  to  fit  a  prior  distribution  to  R,  The  change  of  variable  6  *- 
then  leadR  to  a  prior  distribution  cn  0.  For  that  matter,  since  the  posterior 
distribution  is  eventually  of  interest  the  posterior  distribution  of  6  can  be 
obtained  by  using  the  above  change  of  variable  in  the  posterior  distribution 
of  R.  Another  attributes  situation  arises  when  test  time  T  is  fixed  and  the 
equipment  is  immediately  repaired  when  a  failure  occurs.  Then  (because  of 
the  exponential  assumption)  the  number  of  failures  occurring  in  time  T  is 
Poisson  and  the  K-R  or  Tucker  method  may  be  used  to  fit  g(0). 

The  variables  data  situation  occurs  when  the  actual  failure  times  are 
available.  These  times  occur  either  by  agreeing  to  stop  testing  after  either 
a  fixed  number  of  failures  have  occurred  or  after  a  fixed  time  has  elapsed. 

The  former  case  is  called  a  censored  test  and  the  latter  is  called  a  trun¬ 
cated  ^est.  Often  the  failure  times  themselves  are  not  available  but  the 
pair  (8  ■  observed  MTBF,  K  *  number  of  failures)  is  available.  It  really 
doesn't  matter  for  all  members  of  the  class  of  sufficient  statistics  for  8 
(to  which  6  and  the  failure  times  belong)  result  in  the  same  posterior  dis¬ 
tribution  for  0.  Under  the  exponential  assumption  the  distribution  of  8  is 
gamma  and  the  X-R  method  can  be  used. 

Thus,  virtually  all  forms  of  failure  data  can  be  used  to  fit  prior  dis¬ 
tributions.  However,  in  this  study,  only  the  following  two  types  of  variables 
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data  were  found  to  be  of  practical  int« rest: 

Type  1  -  the  observed  data  is  X,  tne  number  of  failu: :s 
in  a  fixed  time  T. 

A 

Type  2  -  the  observed  data  is  e,  the  sample  mean. 

Both  Type  1  and  Type  2  data  were  used  in  this  study,  with  v.ie  methods 
of  fitting  priors  and  the  results  given  in  Section  4  of  this  re  ort.  The 
discussion  of  quantitative  data  requirements  in  the  next  section  is  re¬ 
stricted  to  just  these  two  types  of  data. 

5.2  THE  AMOUNT  OF  DATA 

In  order  to  establish  quantitative  requirements  for  fittir <  a  prior 
distribution,  the  following  information  is  necessary: 

1)  The  type  of  data  used. 

2)  The  method  used  to  fit  the  prior  distribution. 

3)  The  criteria  established  tc  define  "suitable"  results  for  each 
method  and  type  of  data. 

Because  of  the  dependence  of  the  data  requirements  on  the  above,  each  of 
the  cases  of  Section  4  must  be  analyzed  separately.  For  each  case  (i.e., 
method  and  data  type),  a  sensitivity  analysis  will  be  carried  out,  and  data 
requirements  will  be  set  by  applying  suitability  criteria  to  the  results  ui 
the  analysis. 

5-2.1  THE  AMOUNT  OF  DATA -FAMILY  UNSPECIFIED 

The  only  method  used  in  this  study  for  fitting  the  prior  distribution 
when  the  family  is  unspecified  is  the  Krutchkoff -Rutherford  method.  The 
K-R  method  was  used  on  field  data  in  this  study,  and  the  results  are  sum¬ 
marized  in  Section  4.3.1.  The  purpose  of  this  section  is  to  derive  suitable 
quantitative  data  requirements  for  applying  the  K-R  method. 

A  The  type  of  data  necessary  for  using  the  K-R  method  is  observations 
(bi,K.)  on  n  units.  Since  the  data  is  vf  Type  2,  it  is  necessary  to  set 
minimal  quantitative  requirements  for  the  following: 

1)  the  unit  sample  size,  n. 

2)  the  number  of  failures,  i-1,  . . . ,  n. 

For  the  purpose  of  determining  these  data  requirements,  the  following 
computer  simulation  was  carried  out: 
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n  prior  distribution  ,i(o;  was  assumed  so  be  inverted  gamma  with  param¬ 
eters  a  s-  15J0,  a.  *  5*  ?rom  Section  y.i,  it.  car.  be  seen  that  the  condition 
X  >  4  guarantees  that  the  first  four  moments  of  g(o)  are  finite.  Forty-twc 
pairs  (n,K)  were  selected,  by  combining  6  values  of  n  (n  =  10,20,30,50,100, 
200;  with  each  of  7  valuers  of  K(K  =  5,10, 20, 30, 50,100,200; .  For  each  (n,K) 
pair,  n  random  values  of  d,  each  based  on  K  failures,  were  drawn.  (The 
method  of  drawing  the  random  b's  is  explained  explicity  in  the  next  section.) 
Then  with  the  n  6's,  the  first  four  moments  of  g(b)  were  estimated  as  ex¬ 
plained  in  Section  4.3-1,  and  the  K-R  method  was  applied.  Tne  criterion  for 
"suitable''  data  was  (roughly)  that  an  (n,K)  combination  is  suitable  if  the 
corresponding  simulated  data  leads  to  a  Pearson  curve  reasonably  "close"  to 
the  original  inverted  gamma  prior  distribution  g(y).  One  would  expect  the 
estimated  prior  distribution  to  become  "closer"  to  the  true  prior  as  both  n 
and  K  increase. 


However,  the  result  of  the  simulation  was  that,  in  all  42  cases,  the 
K-R  method  was  not  even  usable.  In  all  rases,  the  estimate  of  the  fourth 
central  moment  was  negative,  thus  making  the  K-R  method  impossible  to  use 
(for  in  all  distributions,  central  moments  of  even  order  are  positive). 

Cince  the  42  cases  are  well  representative  of  the  practical  range  of 
a.ailable  data,  and  the  specific  inverted  gamma  prior  selected  is  very 
"typical,"  the  results  lead  to  the  following  conclusion:  No  data  require¬ 
ments  can  be  set  for  the  K-R  method,  since  the  method  generally  fails  for 
practical  ranges  of  data. 
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In  this  study,  two  cases  of  fitting  a  prior  with  the  family  specified 
were  considered: 

1)  Inverted  gamma  prior,  observed,  random  variable  is  b .  In  this  case, 
minimum  data  requirements  must  be  set  for  (i)  the  sample  size,  n, 
and  ( ii >  the  number  of  failures  ,  i«l,  ...,  n. 

2)  Inverted  gamma  prior,  observed  random  variable  is  X,  the  number  of 
failures  in  time  T.  In  this  case,  a  minimum  data  requirement  must 
be  set  only  for  the  sample  size  n.  (Tne  number  of  failures  in 
this  case  .U  the  random  variable,  not  a  parameter  that  can  be  con- 
trolleu  in  a  test.) 

5. 2.2.1  DATA  REQUIREMENTS  WHEN  THE  OBSERVED  STATISTIC  IS  « 

The  general  framework  is  that  there  is  an  inverted  gamma  prior  g(b)on 
tile  mean  b  of  an  exponential  time-to-failure  distribution.  Suppose  tiiat  for 
each  of  n  equipments,  K  failure  timesAare  observed  and  u  £  is  computed.  One 
then  fits  g(b)  using  the  n  values  of  y.  (This  is  the  Type  2a  data  discussed 
in  Section  4.)  To  test  the  sensitivity  of  the  fit,  one  can  try  fitting  tbe 
sample  data  to  alternate  priors,  and  can  repeat  the  test  for  different  sets 
of  (n,K).  If  g(b)  is  the  true  prior  distribution,  then  intuitively,  as  n 
and  K  are  increased,  the  test  should  accept  g(b)  as  the  prior,  and  reject 
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the  alternate  distributions,  with  increasing  frequency.  In  order  to  get 
concrete  results,  i.e.,  specific  values  of  n  and  K  for  which  discrimination 
between  the  true  and  alternate  priors  is  "good,"  it  was  necessary  to  conduct 
simulations  for  specific  cases. 


The  simulation  program  reads  in  the  parameters  a  and  \  of  g(0),  along 
with  the  parameters  of  the  alternate  distributions  which  were  selected  to  be 
the  Wflbull  and  lognormal  distributions  with  the  same  mean  and  variance  as 
the  inverted  gamma.  The  program  also  reads  in  the  desired  values  of  n  and  K 
(n  =  10,20,30,50,100,200;  K  =  10,20,30,40).  For  each  of  the  24  pairs  (n,K), 
the  following  steps  are  carried  out: 


1)  Draw  n  random  observations  from  g(6). 

2)  For  each  of  the  n  e's,  draw  K  random  failure  times  x^,  ...,  from 

f(x|e)  and  compute  " 

§  -  £  x  /K. 
i=l 

3)  Using  the  n  e's,  perform  tests  to  accept  or  reject  each  of  the 

3  candidate  priors,  at  both  the  *90  and  .95  confidence  levels 

(6  )C  tests  in  all). 


The  whole  experiment  is  repeated  400  times,  and  the  final  output  is  a  table 
th*t  for  e«ch  n  and  each  K,  prints  out  the  number  of  the  4Q0  that  resulted  in 
acceptance,  for  each  of  the  6  X2  tests. 


In  order  to  take  the  simulated  a,  tests,  it  was  necessary  to  compute  the 
percentage  points  of  the  marginal  distributions  corresponding  to  the  inverted 
gamma,  lognormal,  and  Weibull  prior  distributions.  Since  all  three  marginals 
are  analytically  intractable,  an  auxiliary  simulation  was  written  to  do  the 
following: 


1)  10,000  random  samples  aie  drawn  from  g(0). 

2)  From  each  of  these,  h.  random  failure  times  from  an  exponential 
with  mean  6  are  drawn,  and  a  0  is  computed. 

3)  The  10,000  &'s  are  ordered. 

4)  The  100t4,  200^,  . ..,  9900th  s  are  printed  out  in  a  deck  to  be 
used  in  the  main  program  as  an  approximate  "look-up"  table  of  the 
99  percentage  points  of  fK(«). 

2 

The  X  tests  are  carried  out  in  the  main  program  as  follows: 


1)  When  the  n's  are  read  into  the  program,  corresponding  Cq'b, 
denoting  the  desired  number  of  X*  cells,  and  corresponding  .90 
and  .95  level  values  (Cn-1  degrees  of  freedom),  are  also  read  in. 
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2)  For  a  given  (n,K),  class  marks  are  selected  from  the  look-up 
table  so  that  the  expected  number  of  observations  of  q  in  each 
cell  is  n/c^  under  the  hypothesis  that  the  marginal  distribution 

is  the  one  corresponding  to  the  inverted  gamma  prior  g(b). 


3)  Using  the  look-up  tables  corresponding  to  the  K^k-crder  marginals 
when  the  prior  is  Weibull  and  lognormal,  the  expected  number  of 
observations  in  each  cell  under  the  Weibull  and  lognormal  hypotheses 
are  found. 


4)  The  number  of  values  of  0  falling  into  each  of  the  cells  are 
tallied  and 


X 


Cn  (0  -E ,)2 

y  -i— L_ 

k  Ei 


is  computed  for  all  3  cases,  where  i8  the  observed  number 
of  &K's  in  the  itb  cell,  and  the  expected  number  under 

a  given  hypothesip.  In  each  case,  X^  is  compared  to  both 
the  .90  and  .95  X  values,  and  the  hypothesis  is  rejected 
if  X2  is  too  large. 


^  For  one  of  the  400  iterations  of  the  program,  a  complete  output  of  the 
X  tests  was  printed  out  for  each  of  the  24  (n,K)  combinations.  Table 
5. 2. 2. 1.1  shows  the  output  for  the  case  n  *  100,  h  •  bO.  (The  inverted 
gamma  parameters  are  a  *  1500,  X  »=  4).  Note  that  in  this  particular  case 
the  X2  test  is  passed  at  both  levels  under  the  hypotheses  that  the  prior 
is  inverted  gamma  or  lognormal,  and  i6  failed  at  both  levels  under  the 
Weibull  hypothesis. 


Due  to  the  fact  that  the  simulation  of  the  percentage  points  of  the 
marginal  distributions  fK(y)  requires  a  large  amount  of  computer  time,  the 

study  was  restricted  to  one  inverted  gamma  prior  with  parameters  a  *  1500, 

A.  *  4.  The  alternate  distributions  considered  were  the  lognomJL  with 
parameters  u  «  500,  0  •  353*55,  and  the  Weibull  with  parameters  a  ■  8600.7, 
B  *  1*4355*  All  three  distributions  have  the  same  mean  and  standard  devia¬ 
tion,  namely,  u  »  500,  o  ■  353*55*  Plots  of  the  three  density  functions 
arc  3hown  in  Figure  5*2. 2. 1.2. 

The  output  of  the  simulation  is  shown  in  Tables  5*2. 2. 1.3. a  and  g 
5*2. 2.1. 3* b.  For  each  pair  (n,K),  400  simulations  were  run,  taking  2  x 
tests  (.90  and  .95  level)  for  each  of  the  above  3  priors.  The  number  of 
tests  passed  out  of  400  are  given  for  the  .90  level  tests  in  Table 
5. 2. 2. 1.3. a  and  for  the  -95  level  tests  in  Table  5*2.2.1.3*b. 

Beferring  to  the  above  tables,  the  following  rationale  for  selecting 
suitable  (n,K)  pnlre  is  given: 
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TABI£  5. 2. 2. 1.1  EXPECTED  .AND  OBSERVED  VALUES  FOP  X  TESTS 
UNDER  EACH  OF  3  HYPOTHESIZED  PRIORS 
n  =  100 
K  ~  40 


2 

If  N  X  tests  are  run  at  significance  level  p,  then  the  distribution 
of  the  number  of  tests  passed,  under  the  null  hypothesis  Hq  that  the  samples 

come  from  |he  distributior  being  fitted,  is  binomial  with  mean  u  =■  Np  and 
variance  a  ■  Np(l-p). 

2 

In  case  a,  N  *  400,  p  *  .90,  so  that  u  *  36Q,  0  ■  36. 

o 

In  case  b,  N  «*  400,  p  »  .95,  so  that  M.  *  330,  0"  *  19* 

In  each  case,  letting  the  random  variable  X  be  the  number  of  tests  passed, 
the  distribution  of  (X  -u)/o  is  approximately  distributed  N(0,l)(unit  normal). 
If  we  solve  the  equations 


-1-28 


-1.64 


for  x  ^  and  x  QC.,  then 
.yu  *y? 


Pr(x  *  x  ^JHq)  *  .90  and  Pr(x  a  x t95lHp)  -  .95- 

Referring  to  Tables  5.2.2.1.3-a  and  5*2.2.i.3.b,  the  following  criterion 
seems  reasonable:  "Suitable"  pairs  (n,K)  for  fitting  the  prior  are  those 
for  which  the  number  of  "passes"  exceed  x  ^  (n:  x  ^  for  the  test  of  the 

inverted  gamma  prior  (the  true  prior)  and  are  less  than  x  ^  (or  x  ^)  for 

the  tests  of  the  lognormal  and  Weibull  priors  (the  false  hypotheses).  The 
computed  values  of  x  ^  and  x  ^  (rounded  to  the  nearest  integer)  are: 


X.90 

X*95 

Case  a  (.90  level) 

352 

350 

Case  b  (.95  level) 

374 

373 

Looking  at  Table  5 >2. 2. 1.3. a,  one  can  note  the  following: 


1) 


Excluding  the  case  where  n  **  10  (all  K),  and  the  two  cases 

n  ■-  20,  K  »  5  and  n  =  20,  K  «  10,  the  number  of  passes 

fcr  the  lognormal  and  Weibull  tests  ,.s  always  less  than 

b'th  x  q-  and  x  (352  and  350),  whereas,  for  the  in- 
.yj  *yp 


\arted  gamma  test,  only  4  cases  pass  less  than  352  times, 
cnd  only  1  less  than  350  times. 


2)  ?or  each  fixed  K,  the  number  of  times  the  lognormal  and 
Weibull  tests  pass  monotonically  decreases  as  n  increases 
(with  one  exception).  The  same  holds  true  (with  2  exceptions) 
for  fixed  n  f  10  and  increasing  K. 
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TABLE  5. 2.2.1. 3. a.  NUMBER  OF  TESTS  PASSED  (OUT  OF  400) 
UNDER  EACH  OF  3  HYPOTHESIZED  PRIORS,  TAKEN  AT  .90  LEVEL 


10 

1  20 

i 

30 

50 

100 

200 

Inverted  Ganna 

357 

375 

364 

351 

367 

360 

Lognormal 

357 

359 

342 

311 

308 

268 

Weibull 

357 

320 

272 

175 

84 

14 

Inverted  Ganna' 

364 

365 

364 

345 

357 

351 

Lognormal 

364 

351 

347 

! 

302 

251 

210 

Weibull 

364 

299 

, 

233 

131 

31 

0 

Inverted  Gamma 

356 

359 

365 

358 

359 

359 

Lognormal 

356 

330 

319 

313 

239 

to 

w\ 

r-t 

Weibull 

356 

1 

251 

17  2 

86 

9 

0 

Inverted  Gamma 

356 

367 

359 

365 

367 

349 

Lognormal 

356 

318 

298 

261 

189 

85 

Weibull 

356 

244 

161 

75 

7 

0 

m# ! 


An  investigation  of  Table  5*2.2.1.3-b  yields  essentially  the  same  re¬ 
sults  a3  above. 


As  can  be  seen  in  the  tables,  the  discrimination  between  the  true  and 
alternate  priors  increases  rather  dramatically  as  n  and  K  increase.  It  can 
be  seen  why  the  Weibull  distribution  fares  '‘worse"  than  the  lognormal  by 
looking  at  the  plot  of  the  three  densities  (Figure  5. 2. 2. 1.2).  Obviously, 
the  less  the  alternate  prior  "looks  like"  the  true  prior,  the  better  the 
discrimination  between  the  two. 


On  the  basis  of  the  above  results,  the  following  conservative  data 
requirements  are  set  for  fitting  an  inverted  gamma  prior  when  the  observed 

statistic  is  (h  n  >  20  and  K  >  10.  We  must  keep  in  mind,  of  course,  that 
the  requirements  are  based  on  experience  with  one  particular  prior  (inverted 
g&mna  with  «  *  1500,  X  <*  4)  and  two  particular  alternate  priors.  It  would 
be  of  interest  to  repeat  the  simulation  for  other  cases.  It  is  felt,  however 
that  the  case  considered  is  fairly  representative  and  that  the  results  ob¬ 
tained  from  it  are  generally  applicable. 


The  only  "problem"  with  the  analysis  described  qbove  is  that  it  is 
limited  to  cases  where  K  is  fixed  in  each  sample  of  b'e  [Type  2a  lata). 
When  the  simulations  were  first  undertaken,  it  was  thought  coat  if  a  pair 


(n,K)  was  suitable  for  a  prior  fit,  then  data  of  the  form  ,Ki),i»l,  .  ,.,n, 

all  K.  is  k}  would  also  be  suitable.  This  is  not  the  case,  however,  since 
when  the  K's  are  different  (Type  2b  data),  one  must  use  the  mixed  population 
model  for  the  marginal  distribution.  As  explained  in  Section  4.3.2,  this 
leads,  in  actual  practice,  to  a  problem  concerning  the  degrees  of  freedom 
for  the  X2  test,  due  to  n  being  too  small  and/or  having  too  much  "variety" 
in  the  K. 's.  In  the  above  simulation,  this  problem  was  avoided  by  having 
K's  for  each  sample.  The  basic  quantitative  data  requirement  when  the  ' s 

are  different  is  that  derived  in  Section  4.3.2:  c  a  t  +  4,  where  c  is  the^ 

2 

number  of  X  cells,  and  t  is  the  number  of  distinct  K^/s.  Since  we  require 


n 

5 


2 


we  can  rewrite  the  requirement  as  t  +  s 


This  requirement  was  not  achieved  in  this  report  as  can  be  seen  in 
Section  4.3.2.  Two  renarks  are  in  order.  First,  no  matter  how  many  distinct 
K^s  there  are  (even  if  t-n)  if  the  ^  are  all  large,  i.e.,  greater  than 
or  eqqal  to  20,  then  the  prior  distribution  can  be  fitted  directly  to  the 
observed  f's.  Secondly,  when  n  is  relatively  large  with  respect  to  t  the 
mixed  model  can  be  fitted.  There  will  be  cases  for  large  n  when  t  is  rela¬ 
tively  small. 

5. 2. 2. 2  DATA  REQUIREMENTS  WHEN  THE  OBSERVED  STATISTIC  IS  THE  NUMBER  OF 
FAILURES  IN  TIME  T 


Sections  4. 2. 2.1  and  4. 2. 3-2  discuss  the  method,  and  Section  4.3,2  the 
results,  of  fitting  an  inverted  gamma  prior  distribution  when  the  observed 
statistic  is  the  number  of  failures  X  in  a  fixed  time  T.  The  purpose  of 


TABLE  5. 2. 2.1.3.  b.  NUMBER  OF  X  TESTS  PASSED  (OUT  OF  400) 
UNDER  EACH  OF  3  HYPOTHESIZED  PRIORS,  TAKEN  AT  .95  LEVEL. 


X 

i 

10 

20 

'  ' 

30 

50 

100 

200 

Inverted  Gamma 

394 

390 

382 

379 

381 

376 

✓ 

Lognormal 

394 

377 

373 

343 

348 

299 

Weibull 

394 

349 

318 

233 

124 

22 

Inverted  Gasnna 

393 

385 

381 

381 

376 

377 

10 

Lognormal 

381 

373 

372 

343 

306 

259 

Weibull 

39 3 

335 

286 

186 
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1 

Inverted  Gamma 

392 

379 

380 

378 

382 

379 

20 

Lognormal 

377 

364 

357 

343 

290 

203 

Weibull 

392 

286 

229 

127 

17 

0 

Inverted  Gamma 

388 

380 

376 

383 

387 

364 

40 

Lognormal 

388 

360 

338 

320 

250 

130 

Weibull 

376 

284 

217 

114 

17 

0 

this  section  is  to  find  the  minimum  unit  sample  size  n  required  to  get  a 
"suitable"  prior  fit.  The  analysis  to  determine  the  minimum  n  was  based  on 
the  sensitivity  study  described  belov,  which  is  analogous  to  the  one  de¬ 
scribed  in  the  previous  section.  The  restriction  on  T  Is  that  it  be  large 
enough  so  that  the  number  of  failures  is  variable  enough  to  take  the  X2  test. 

The  inverted  gamma  prior  parameters  selected  for  the  simulation  were 
or  *  1*000,  X  »  3*  As  in  Section  5*2. 2.1,  two  alternate  priors  were  selected 
having  the  same  mean  and  variance  as  the  above  prior:  a  Weibull  distribution 
with  parameters  a  ■  2000,  8  ■  1  and  a  lognormal  distribution  with  parameters 
Vi  ■  2000,  o  •  2000.  The  fixed  time  T  was  chosen  to  be  T  *  1*000  hours.  Six 
values  for  the  sample  size  n  were  studied:  n  »  10,20,30,50,100,200.  For 
each  n,  n  random  values  of  X,  {x^,  ...,  xnj  were  drawn,  and  x2  tests  taken 

against  the  (true)  Inverted  gamma  prior  and  the  alternate  hypothesized 
priors.  Each  random  value  of  X  was  drawn  by  taking  a  random  value  of  0 
from  the  inverted  gamma  prior  g(8),  and  then  taking  random  failure  times 
t^,tg, ....  from  the  exponential  distribution  with  mean  6 .  The  value  of  X 

K 

is  then  taken  to  be  the  largest  value  of  K  for  which  T'  t.  <  T. 

Pi  ^ 


Table  5. 2. 2. 2.1  shows  the  results  of  one  run  of  the  simulation  for 
the  case  n  ■  200.  The  first  column  of  the  upper  table  lists  the  value,  x, 
that  the  random  variable  X  can  take  on,  and  the  second  column  gives  the  ob¬ 
served  number  of  times  (out  of  200)  that  the  events  X*x  occurred.  The 
last  3  columns  give  the  expected  number  of  times  (out  of  200)  of  the  occur¬ 
rence  of  the  events  X»x  under  the  three  hypothesized  prior  distributions. 
These  expectations  are  easy  to  compute  fear  the  case  of  the  inverted  gamma 
prior,  since  the  marginal  distribution  of  X  can  be  computed  explicitly  by 
Formula  4. 2. 2. 1.2.  For  the  lognormal  and  Weibull  cases,  however,  the  mar¬ 
ginal  distributions  had  to  be  approximated  by  simulation.  In  both  cases, 
10,000  random  values  of  X  were  drawn  (by  taking  random  0's  from  g(0)  and 
obtaining  a  random  value  of  X  as  described  in  the  previous  paragraph) .  Then 
the  relative  frequencies  (out  of  10,000)  of  the  events  X  ■  0,1,2,...,  were 
used  to  approximate  the  marginal  f(x). 

2 

In  the  lower  table  of  Table  5<2.2.2.1,  the  results  of  the  X  tests  are 
given.  Six  degrees  of  freedom  were  used  because  the  data  in  the  upper  table 
was  divided  into  7  X2  ceils.  For  this  particular  case,  the  inverted  gamma 
parlor  passes  the  X2  test  at  both  the  .90  and  .95  level,  while  the  other  two 
hypothesised  priors  fall  both  tests. 


As  in  Section  5*2. 2.1,  the  above  simulation  was  repeated  1*00  times  for 
each  n.  Table  5-2. 2.2.2  gives  the  number  of  times  out  of  400  that  each  test 
was  passed.  The  criterion  for  "suitable"  values  of  n  is  the  same  as  in 
Section  5. 2. 2.1,  leading  to  the  same  critical  values  x  ^  and  x  ^  as  given 

In  that  section.  As  can  be  seen  in  the  table,  when  n  is  t  30,  the  number 
of  times  the  Inverted  gamma  prior  passes  the  x  test  exceeds  the  critical 
value  350(373)  the  <90(.95)  level,  whereas,  tbs  number  of  passes  for  the 
lognormal  and  Wtlbull  priors  fall  short  of  the  critical  value. 
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Hence,  in  the  case  of  fitting  a  prior  specified  as  inverted  g»unmw.  when 
the  observed  data  is  the  number  of  failures  in  time  T,  the  following  conser¬ 
vative  data  requirement  is  set:  the  sample  size  n  must  be  at  least  30. 


TABLE  5. 2. 2. 2.1  EXPECTED  AND  OBSERVED  VALUES  FOR  X  TESTS 
UNDER  EACH  OF  3  HYPOTHESIZED  FRIORS  (n=200) 


Value  x 
of  X 

Observed  No. 

0*  lines  X  -  x 

Expected  No. 
for  Inverted 

ustma  Prior 

Expected  No. 
for  Log  Normal 
Prior 

Expected  No. 
for  Weibull 
Prior 

0 

33 

25.00 

25.16 

28.76 

1 

43 

37.50 

35.06 

34.20 

2 

41 

37.50 

32.08 

26.80 

3 

29 

31.25 

24.92 

20.68 

4 

i3 

23.41 

19.60 

14.52 

3 

16 

16.41 

15.80 

10.20 

6 

8 

10.94 

11.34 

7.70 

7 

11 

7.03 

8.24 

7.12 

Over  7 

6 

10.94 

22.06 

44.62 

Computed  )T 
.90  Level 
.95  Lml 


Inverted 


8.964 

10,645 

12.592 


Lognormal 

Weibull 

44.196 

20.034 

10.645 

10.645 

12.592 

12.592 

TAJEJIE  5*2. 2. 2. 2  NUMBER  OF  X2  TESTS  (OUT  OF  *00) 

UNEER  EACH  OF  3  HYPOTHESIZED  ERICKS 


.90  Inverted  Genoa  352  342  353  355  353  367 

Level  Lognormal  352  332  328  260  183  57 


Weibull  373  332  263  46  11  0 

.95  Inverted  Guam  390  3«0  383  373  376  387 

Level  Lognormal  373  364  352  307  228  8? 

Weibull  373  364  304  80  16  0 


SECTION  6.0  ANALYSIS  FOR  DATA  COMBINATION 


In  fitting  prior  distributions  in  this  study  different  classes  of 
identical  equipments  were  available.  In  fact,  seventeen  data  sets  were 
used.  It  turns  out  that  some  care  must  be  exercised  in  defining  what  is 
meant  by  a  prior  distribution  on  MTEF,  say  g(6).  All  p.d.f.'s  must  exhibit 
some  form  of  homogeneity.  That  is,  limitations  are  placed  on  the  "reasons" 
why  the  variable  of  interest  varies.  These  reasons  are  commonly  called 
assignable  causes  of  variation.  As  a  beginning,  g(B)  has  been  restricted 
as  foli 

The  random  MTBF'6  (e's)  belonging  to  a  particular  g(e) 
must  be  6's  on  a  given  type  equipment  built  to  the  same 
design  specifications  by  a  particular  manufacturer. 

This  is  somewhat  restrictive  because  it  apparently  places  equipment 
of  the  same  design  but  built  by  different  manufacturers  in  different  g(9)'s. 
It  also  places  similar  equipment,  say  computers,  but  witn  different  designs, 
e.g. ,  different  memory  sire,  in  different  g(e)'s.  One  might  ask,  couldn't 
some  of  this  data,  say  different  computer  types,  be  combined  into  one 
prior  distribution?  The  answer  is  yes,  but  with  this  qualification:  every 
assignable  cause  of  variation,  e.g.,  different  manufacturers,  is  a  piece  of 
prior  information  and,  if  possible,  should  be  exploited.  Combining  data 
into  one  prior  which  have  assignable  and  identifiable  causes  for  having 
different  MPBF’s,  in  general,  increases  the  variation  in  the  prior  distribu¬ 
tion.  On  the  other  hand,  fitting  a  large  number  of  prior  distributions  is 
a  costly  process  and  it  is  worthwhile  to  be  able  to  "relate"  prior  distribu¬ 
tions  on  similar  equipments  even  though  they,  the  prior  distributions,  are 
not  combined.  For  example,  two  computers  of  similar  design,  say  different 
memory  sizes,  might  ha,fe  prior  distributions  which  are  relatable,  though  it 
might  not  he  wise  to  combine  them.  Thus,  when  the  prior  distribution  of 
the  one  is  fitted,  the  prior  distribution  of  the  other  is  known.  For 
example*  let  these  computers  be  called  c^,  cg  with  prior  distributions  g^fc) 

and  82(e)-  It  might  turn  out  that  g^($)  and  82(e)  belong  to  the  same  family, 

«.g. ,  Inverted  gamma,  and  hence  may  Le  related  by  some  transformation.  If 
this  transformation  is  known  then  having  fitted  81(e)(g2(b))# ^(wHg-^fe))  i8 

known.  In  Section  9.3  of  the  Appendix,  a  particularly  simple  transformation 
was  considered,  namely, 

e2  -  K6X  K  >  0. 

Then,  if  is  inverted  gamma  with  parameters  (ar^,  X^),  *2  is  inverted  gamma 

»itb  parameters  (og  *  «^K,  lg  *  X^).  Thus,  knowing  (or^,  X^)  and  the  constant 
K,  the  parameters  of  the  second  prior  distribution  are  known. 
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One  particularly  goad  idea  for  estimating  K,  which  could  not  be  verified 
oecause  of  lack  of  data  on  two  similar  equipments,  is  that  K  might  be  well 
represented  by 


K  =  fa  M 
P2/  PX 


where  ep  is  the  predicted  MTBF  on  the  1th  equipment.  This  idea  should  be 


investigated  in  the  next  study  phase. 


SECTION  7-0  ROBUSTNESS  ANALYSIS 


SUMMARY 


Objective 

The  objective  of  this  section  was  to  investigate  the  effects  of  errors 
in  estimating  the  scale  and  shape  parameters  of  inverted  gamma  prior  distri¬ 
bution  on  the  posterior  inverted  gamma  distribution.  The  effects  were 
measured  in  terms  of  changes  in  the  mean,  5^  lG^k,  90  and  95ttl  percentiles 
of  the  posterior  inverted  g*anw.  distribution.  For  measuring  effects  of 
errors  in  estimating  the  scale  parameter  a,  a  shape  parameter  '.*=3  was 
assumed.  For  measuring  effects  of  errors  in  estimating  the  shape  parameter 
^  scale  parameter  o*>200  was  assumed.  The  posterior  distribution  depends  on 

6  and  three  values  were  chosen:  &-50,  100,  200.  These  choices  were  reason¬ 
able  in  terms  of  the  value  of  o,X  assumed;  K  values  chosen  were  5,10,20,30. 

Results 


The  effect  of  errors  in  estimating  the  scale  parameter  of  the  prior 
Inverted  gamma  distribution  were  practically  negligible  on  the  selected 
percentile^  of  the  posterior  inverted  ganaa  distribution  for  all  three 
values  of  9  and  for  K  *  20.  The  effects  were  more  noticeable  for  K*10  and 
quite  pronounced  for  K»5-  This  ties  in  with  the  results  of  Section  5-0 
regarding  data  requirements  although  arrived  at  in  a  different  way. 

The  effects  of  errors  in  estimating  the  shape  parameter  1  in  the  prior 
inverted  ganaa  distribution  were,  as  in  the  above  case,  practically  negligi¬ 
ble  on  the  selected  percentiles  for  K  *  20.  For  K  »  5,10  the  effects  were 
quite  pronounced.  This  result  again  ties  in  with  the  results  of  Section  5.0 


It  had  originally  been  intended  to  select  several  Bayes  reliability 
demonstration  test  methods  and  investigate  their  robustness  with  respect  to 
errors  in  estimating  the  parameters  and  family  of  the  prior  distribution. 
However,  only  two  methods  (Bibliography  #2,  #59)  are  available  at  this  time 
and  neither  has  yet  gained  anything  near  acceptance  in  the  field  of  reliability 
testing.  For  these  reasons,  it  was  decided  to  abandon  this  approach  and  take 
another  tack  which,  hopefully,  would  be  of  mere  use  to  the  reader. 

The  ultimate  use  of  the  prior  distribution  is  to  supplement  it  with  ob¬ 
served  data  so  that  a  posterior  distribution  may  be  calculated.  In  Bhort, 

Bayes  reliability  tests  will  use  the  posterior  distribution  in  one  forr.  or 
another.  Thus,  it  was  decided  to  investigate  the  sensitivity  of  posterior 
distribution  to  errors  in  estimating  the  parameters  of  the  prior  inverted 
gamma  distribution;  there  was  not  time  to  investigate  differences  in  families. 


1  r  (u,  X)  are  the  scale  and  shape  parameters  of  a  prior  inverted  gaxma 
distribution,  then  the  posterior  distribution  has  parameters  {«  KB,  X  •*  K) 
with  mean 

e(*!§)  -  jSgf  (7.1) 

!iow,  if  (a, X)  were  estimated  incorrectly  as  say,  a1  “  or  ♦  X*  ■  X  ♦  Al  then 

(7-z> 

and  this  ^ould  be  the  posterior  mean  used.  .However,  for  large  K  both  E(ei$) 
and  E*(t>lo)  approach  the  same  limit,  i.e.,  B,  and  do  not  differ  by  much.  The 
percentiles  of  the  posterior  distribution  are  not  available  in  closed  form 
but  can  easily  be  obtained  by  computer.  Figures  7.1  through  7.12  present 
the  behavior  of  the  posterior  Inverted  gamma  distribution  (in  terms  of  its 
mean,  5bb,  10th,  90bh  and  95th  percentiles)  for  a  fixed  value  of  X  •  3  as  a 
function  of  the  scale-  parameter  a.  The  X  ■  3  vas  selected  because  it  repre¬ 
sented  a  reasonable-  magnitude  for  the  nhape  parameter.  Since  the  posterior 

distribution  depends  on  6,  three  &'a  were  selected  for  each  K;  50,  100,  200. 
The  K’s  selected  were  5,  10,  20,  30-  Figure  7.x  shows  that  the  mean  add 
percentiles  of  the  posterior  distribution  are  rather  insensitive  to  a  (i.e., 
the  lines  are  almost  vertical)  for  X  *  30,  §  «  200.  Hicb  the  same  is  true 
for  all  K  *  30,  i.e.,  for  Figures  7.1,  7-2,  7>3-  When  K  -  10  (Figures  7-7, 
7-6,  7-9)  the  mean  and  percentiles  become  quite  sensitive  to  differences  in  u. 
Note  that  the  effects  of  b  are  in  the  translation  of  the  mean  and  percentiles 
along  the  B  axis  and  not  in  the  slope  of  the  lines.  Thus,  far  X  ■  3  the 
errors  of  estimating  a  incorrectly  do  not  appear  serious  until  K  gets  small. 
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Comparing  Figures  7-1  through  7 *12  with  Figure  7-0  shows  that  though  the 
prior  distribution  i6  highly  sensitive  to  differences  in  a  (for  X  =  3)  the 
posterior  distributions  do  not  exhibit  this  characteristic.  The  posterior 
distributions  are ,  of  course,  sensitive  to  differences  in  a  +  but 

Kfe  is  always  known  and  usually  large  compared  to  u  so  that  errors  in  esti- 
mating  a  are  not  too  serious. 


A  study  similar  to  the  previous  discussion  was  done  for  the  same 

&  *  5Q,  100,  200,  K  >*  5 ,  10  ,  20  ,  30  but  this  time  a  a  200  was  held  fixed  and 
1  varied.  Figure  7*13  shows  the  prior  distribution  mean  and  percentiles. 
Since  the  shape  parameter  of  tae  posterior  distribution  (Figures  7-1^  through 
7-25)  is  (X  +  K),  it-  is  no  surprise  that  the  sensitivity  of  the  mean  and  per¬ 
centiles  to  errors  in  X  is  greater  for  smaller  K.  For  large  K,  (X  +  k)  is 
dominated  by  K  for  the  values  of  X  expected  to  occur  in  reliability.  Put 
another  way,  suppose  X  »  5  in  the  prior  distribution  but  that  an  0O70  error 
was  made  and  X1  «  8  was  estimated.  For  K  =  3 0  then  the  true  shape  parameter 
in  the  posterior  distribution  would  be  (5  +  30)  =  35  and  the  one  estimated 
would  be  (8  +  30)  -  38  a  relatively  small  percentage  error  compared  to  the 
orlcunal  error. 
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SECTION  8.0  CGNCUJSIGNS  AND  RECOMMENDATIONS 


8.1  CGNCUJSIONS 

As  a  result  of  this  study  certain  conclusions  appear  inescapable.  We 
list  the  important  ones  now.  First,  fitting  prior  distributions  to  6  *  MTBF 
is  entirely  feasible.  The  techniques  developed  in  this  report;  Indicate  that 
the  mechanics  of  the  fitting  is  a  relatively  inexpensive  process.  Secondly, 
the  amount  of  data  available,  i.e.,  in  existence  today,  which  leads  to  good 
prior  distribution  fits  is  somewhat  limited.  Thus,  the  costs  of  fitting 
prioi-  distributions,  while  not  exorbitant,  will  be  primarily  incurred  in 
acquiring,  either  by  intensive  search  or  by  designed  test,  suitable  data. 
Thirdly,  the  inverted  gemma  prior  distribution,  which  is  mathematically 
quite  tractable,  also  appears  to  graduate  data  quite  nicely  since  In  the 
eight  (8)  situations  where  a  fit  could  be  obtained,  seven  of  them  were 
well  described  by  an  inverted  gamma  prior  distribution. 

In  susmary,  the  development  of  Bayes  reliability  demonstration  tests 
appears  to  be  quite  feasible  from,  the  standpoint  that  the  prior  distribution 
needed  can  be  fitted. 

8.2  RECOMMENDATIONS 

The  following  re  commendations  are  intimately  connected  with  the  conclu¬ 
sions  of  this  report.  First,  it  is  recommended  that  additional  prior  distri¬ 
butions  be  fitted.  That  is,  prior  distributions  should  be  fitted  to  equipment 
different  than  the  type  studied  in  this  report.  This  will  further  test  the 
suitability  of  the  inverted  gamma  prior  distribution.  Secondly,  in  view  of 
the  costs  of  fitting  prior  distribution,  two  areas  should  be  studied. 

i)  Bayes  methods  of  reliability  demonstration  which  do  not  require 
a  prior  distribution  should  be  studied.  Such  a  method  is 
empirical  Bayes.  It  is  doubtful  this  area  will  be  entirely 
fruitful  but  is  certainly  worth  a  look. 

ii)  Methods  of  relating  prior  distributions  on  different,  but 
perhaps  similar,  equipments  so  that  prior  distributions  may 
be  derived  from  one  another  instead  of  fitting  new  prior 
distributions  to  each  one.  An  example  would  be  relating  two 
equipments  by  some  fixed  ratio  of  their  predicted  MTBF's. 

Finally,  the  receptiveness  of  Government  and  Industry  to 
Bayes  Demonstration  tests  should  be  studied  for  even  if 
the  methods  are  developed  and  feasible  the  Bayes  plans  will 
not  return  such  cn  investment  if  they  are  not  used. 

In  summary  then,  it  is  recommended  that  the  next  logical  steps  be  taken 
in  the  development  of  Bayes  reliability  demonstration  tests. 
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SECTION  9*0  AFFENDIX 


9-1  THE  DIVERTED  GAMMA  KttOR  DISTRIBUTION 

At  the  time  of  this  writing  (late  1969)  It  seems  Impossible  to  exaggerate 
the  importance  of  the  inverted  gamma  distribution  in  reliability  estimation/ 
demonstration;  particularly,  when  the  measure  of  reliability  is  6  *  MTBF. 

The  reasons  far  this  importance  follow. 

First,  the  inverted  gsmma  distribution  is  a  two  parameter  distribution 
(it  can  be  made  three  parameter)  and  is  very  flexible.  That  is,  the  inverted 
gamma  distribution  can  be  used  to  graduate  a  wealth  of  empirical  data.  As 
evidence  of  this,  we  offer  the  fact  that  seven  of  the  eight  data  sets  analyzed 
resulted  in  good  inverted  gmma  fits. 

Secondly,  and,  perhaps  more  importantly,  the  marginal  distribution  of 
X  (X  being  either  time  to  failure  or  number  of  failures  per  fixed  time  T)  is 
also  available  in  closed  form.  This  makes  problems  of  fitting  the  prior 
distribution  much  more  tractable. 

Finally,  the  inverted  gaiana  distribution  is  the  natural  conjugate  for  the 
Poisson  process  (l.e.,  for  the  exponential  and  Poisson  conditional  distribu¬ 
tions)  ,  This  means  that  the  posterior  distribution  is  of  the  same  form  as  the 
prior  distribution  (inverted  gamma)  and  hence,  the  posterior  distribution  is 
available  in  closed  farm.  This  makes  computations  for  demonstration  tests 
and  other  analyses  much  more  tractable. 

The  gamma  distribution  in  the  random  variable  X  is  given  by 

g(x)  -  J*  e'0*  «,\,x  >  0.  (9.1.1)* 

The  change  of  variable  6  *  l/x  leads  to 

«(6)  -  e^0  (9-1.2) 

Thus,  the  gamma  and  inverted  gamma  are  related  by  a  reciprocal  transfor¬ 
mation.  This  is  particularly  convenient  since  for  the  Poisson  process  the 
KTBF  is  the  reciprocal  of  the  failure  rate.  Hence,  a  gamma  prior  distribution 
an  the  failure  rate  implies  an  inverted  gamma  prior  distribution  on  the  MTBF 
and  conversely. 


# 

r(A)  is  notation  for  the  gamma  function:  r(X)  *  f  yx“x  e"y  dy. 

0 
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The  characteristic  function  of  the  Inverted  tfunaa  is  useless  fcr  finding 
the  nonents  hut  they  nay  he  found  directly.  The  ncnent  is  given  hy 


m 

E(eK)  -  y*  eKg(e)de 


.  «x 
rrn 


w 

y  «-«/•« 


(9-1-3) 


The  change  of  variable  y  «  a/ 8  leads  to 


E(eK)  -  /“  y^*^"1  e“ydy 


(9*1.4) 


For  (l-K)  >  0  the  integral  in  (9*1*4)  converges  and  because  of  the 
recurrence  r(v)v  •  r(wl) 


E(8K) 


(9*1*5) 


a  U-i) 

i*i 


For  (X-K)  s  0,  i.e.,  X  a*  K,  the  integral  in  (9*1*4)  is  infinite  and  E(8K) 
does  not  exist. 


If  [X]  *  the  largest  integer  was Her  than  X  (e.g. ,  [4*3]  ■  4,  [4.0]  *  3) 
then  all  nonents  up  to  and  including  exist  and  no  nonents  beyond  [X] 
exist.  Using  (9*1*5): 


E(8)  -  « 

r-r 


X  >  1 


E(82)' 


q2 

(X-l)(X-2) 


\  >  2 


(9*1*6) 


(X-l)Z(X-2) 


EXAMPLE:  Suppose,  given  an  inverted 


distribution  vlth  or  *  100  X  *  2.J. 


Then'  £x]  ■  2  and  only  the  first  tvo  nonents  exist.  They  are 


E(e)  -  ™  -  66.7 
1*5 


? .  1 A  AAA 


!  $ ■ 


■ 1 


13,333*3* 


It  Is  of  no  consequence  to  this  study  that  the  moments  do  not  all  exist. 
Other  than  this  the  inverted  gamma  is  well-behaved  and  even  if  the  mean  does 
not  exist,  the  mode,  median  and  all  quantiles  do.  The  mode  is  given  by 


Mode 


ot 


(9.1.7: 


9.2  THE  MARGINAL  DISTRIBUTIONS 

In  this  section  the  marginal  distributions  (for  an  inverted  gamma  prior 
distribution)  will  be  derived  for  two  conditional  distributions. 

1)  The  ?ol860a  distribution  when  number  of  failures  occurring 
in  fixed  time  T  is  the  observed  random  variable. 


11)  The  cnimn  distribution  when  the  time  to  failure  is  exponential 
and  hence,  the  observed  random  variable  is  sample  KTBF  tJ. 

For  case1  i)  above 


f(x)  f(x|e)g(8)de 


a 

T(x) 


9  ’{X+1)  e  -a/e  I  d0 


] 


.  «  T*  J 


X 
r 

nxy 


ix+x+r 


do. 


With  the  change  of  variable  Z  »  (T*a)/e 

fW  ‘  £  nu 


e“  Z  ZX'*-X  dZ. 


Tbs  Integral  tens  is  rCX**)  so  thst 


f(*) 


RxTiT  (ife)  [fh) 


(92.1) 
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Setting  p  ■  and  q.  ■  ,T  it  is  noted  that  p*q  »  1  and 
THa  Tea 

f*x>  *  ^  ^  x'0’1 . . 

is  the  well-known  negative  binomial  distribution  with 

Mean  *  E(x)  «*  XT 
TT 

Variance  »  *  XT(T+g) 

X  -7"^ 


It  should  be  noted  that  if  we  had  started  with  a  gamma  distribution  an 
failure  rate  and  a  Poisson  conditional  with  mean:  failure  rate  times  T 
Instead  of  t/6  the  same  negative  binomial  distribution  as  (9*2.1)  would 
have  resulted. 


In  case  11 )  above,  it  is  easy  to  show  that  if  time  to  failure  is 
exponential,  then  the  conditional  distribution  of  the  MTBF,  t, 
given  the  MTSF,  6,  is  gasna,  i.e*. 


K 


AVI 

1  (e)**1  e“T" 

m 


here;  X  is  the  observed  number  of  failures  and  8  is  the  true  but  unknown 
MTBF*  Hence,  the  marginal  distribution  of  f  is 


fx(6)  * 


i^e(e)  g($)dd 


:•  .  a 


The  change  of  variable  t  *  lead#  to 

$ 


,<*> .  (  a.V  (AT1  ( K  A 

‘  Ip/  Ug/  Ip/ 


>  0 


which  Is  the  not  so  well-known  inverted  Beta  distribution* 
Jn  any  event 


(9-2*2) 


E(S)  -  £ 

X  >  1 

X  >  2. 

9-3  TRAKSFCRMATIOSIS  OH  THE  INVERTED  GAMMA  DISTRIBUTION 

In  Section  6.0  (Data  Combination)  a  certain  method  of  relating  prior 
distributions  on  "different"  equipments  was  discussed.  The  idea  depends 
on  a  property  of  the  inverted  gamma  distribution  (many  other  distributions 
have  this  same  property)  which  will  be  shown  here: 

If 

g(e)  -  e‘(ui)  e-*^,  (9.3.1) 

ruj 

l.e.,  if  0  is  inverted  gamma  with  scale  parameter  a  and  shape  parameter  X 
then  the  random  variable  y  *  86,  0  <  0  <cc*  is  again  Inverted  gamma  with 
scale  parameter  (ar8)  and  shape  parameter  X.  That  is, 

g(jr)  -  e-“8A  (9.3.2) 

This  result  follows  by  substituting  6  *  /  in  (9- 3* l)  and  multiplying  by 
the  differential  element  dd  *»  dy/fl .  0 

A  change  of  the  form  8  ■  -  c,  0  <  c  has  the  effect  of  introducing 

a  guarantee  time  8c  and  changing  the  scale  parameter  to  ag. 

A  change  of  the  form  0  »  y^  means  that  y  is  not  of  the  Inverted  gamma 
family. 


9.4  IDElfTIFIABILITY  CP  THE  PRIOR  DISTRIBUTION 

The  central  problem  in  estimating  the  prior  distribution  of  MfHF,  say 
g(d),  is  that  random  samples  (to  construct  an  estimate  of  g(8))  are  not 
available  from  g(fi).  Thus,  the  marginal  distribution 

m 

f(x)-J’  f(x|e)g(8)de  (9-4  1) 

0 

must  be  used  to  estimate  g(d)>  Suppose  now  that  the  equation 
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•’flPttffVJIM-*?*.'.-  ,  |T‘-"Dimri  r  r-  .11 '  (rUTrifUtBIH".'  11  U'HIII-  'IS  t 


■■an.  »  ■*« W»>»i  TW 


f(x) 


f(x  le)  gx  (e)  a© 


(9.4.2) 


Implies  that  g^(@)  *  -.■.{&)  •  In  this  situation  f(x)  is  called  identifiable 

with  respect  to  f(xl6).  Identifiability  is  u  crucial  property  far  this  study, 
since  if  (9.4.2)  could  hold  and  g^(&)  t  82(6)*  then  two  different  prior  dis¬ 
tributions  could  lead  to  the  same  marginal  distribution  and  since  tte  Margins! 
distribution  Is  used  to  male  inferences  about  g(&)  it  would  be  impossible  to 
tell  which  (in  the  case  of  non  identifiability)  g(8)  obtained. 


We  will  show  the  identifiability  of  the  marginal  distribution  f(x)  with 
respect  to  three  conditional  densities  (the  Poisson,  exponential,  and  go— r ) 


i)  f(x|«)  »  £  /  foA)  .  ThiB  is  the  case  where  the  number  of  Adheres 

X  • 

occurring  in  fixed  time  T  is  the  observed  random  variable  and  the 
operating  process  is  the  Poisson  process. 


ii) 

iii) 


f(x|e)  -  1/0  e“x/6.  Sere  the  random  variable  X  is  tine  to  failure. 

A 

£ 

e  .  This  is  the  gam*  distribution  of 

sample  MTBF  when  times  tc  failure  are  exponential. 


A  result  of  Teicher  (Reference  7)  is  used  to  show  that  f(x)  is  identifi¬ 
able  w.r.t.  i)  above,  i.e.,  that  different  prior  distributions  cannot  lead  to 
the  same  marginal  distribution.  The  result  is  that  a  family  of  densities 
which  Is  addltively  closed  (a.c.)  is  identifiable.  Additivtly  closed  neens 

f(x|a)  *  f(xj»)  -  f(x|a  ♦  g)  (9-4.3) 

The  *  denotes  convolution.  Sioce  it  is  easy  to  show  that  the  sum  of  two 
Pwisson  variates  is  again  a  Poisson  variate  with  parameter  the  sun  of  tbs 
two  parameters  f(x)  is  Identifiable  in  case  i). 

For  case  ii),  ve  note  that 


f(x)  •  j  l/fi  e “X/eg(e)d9 
o 


(9*4.4) 
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Making  the  change  of  variable  \  =  l/e 

CO 

~  J"  (9. 4. 5) 

o 

Thus,  f(x)  divided  by  X  is  the  Laplace  transform  of  g(x).  If  g(X)  is  a  class 
of  continuous  densities  (e.g.,  the  gamma  family)  then  g(G)is  a  continuous 
class  (e.g.,  inverted  gamma)  and  by  the  uniqueness  theorem  for  Laplace 
transforms  f(x)  is  identifiable. 

Finally,  for  case  iii)  we  again  use  a  result  of  Teicher  (Reference  7)* 
That  is,  for  fixed  K  iii)  is  a  scale  parameter  family  generated  by  f(x|n,l) 
and  a  uniqueness  theorem  for  Fourier  transforms  gives  the  desired  result. 


9*5  LIMITING  BEHAVIOR  OF  THE  MARGINAL  DENSITY  f„(§) 

K 

When  0'  is  the  sample  MTBF  based  on  K  failures,  then  under  the  conditions 
glv«*  belcw  the  marginal  distribution  f^(y)  converges  to  the  prior  distribu¬ 
tion  g(0)  as  K  -»oo  .  This  is  shown  in  the  result  and  proof  9* 5*1  below. 
Figure  9*5*1  shows  a  particular  example:  the  prior  density  g(0)  is  plotted 
along  with  the  corresponding  inverted  Beta  marginal  densities  fK($)  for 

K  =  1,2,3,1,5,10,20,50.  Clearly,  when  K  is  small,  it  is  erroneous  to  fit 
sample  values  of  $  directly  to  g(e).  When  K  is  large,  however,  (any  K  ^  20), 
the  error  is  not  too  bad,  and  sample  based  on  large  K's  can  be  fitted 
directly  to  g(0).  Thus,  when  K  is  "large"  enough,  a  fairly  accurate  method 
exists  for  fitting  the  prior  which  avoids  the  problems  involved  in  using  a 
mixed  population  model  for  the  marginal  distribution. 


9*  5*1  A  RESULT  AND  PROOF 


Consider  a  density  g(0)  defined  on  (0,^9  such  that  sup  g(0)  is  finite, 
and  consider  a  sequence  of  R.V.’s  £t„3  with  densities  given  by 

IV 

CO 

yt)  =  ^g(0)f^.(t 1 0)d9,  K  =  1*2,  ...  where  the 
o 

sequence  {f^}  has  the  property  that  V  fixed  tQ  >  0, 


lim 
K  -*co 


fK(t|0)dt  = 


if  0  <  t 

o 

if  0  >  t 

o 
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LThft  above  condition  is  ^tlefied  in  the  usual  Bayes  *t-up,  vhere  T„  is  (say) 
A  *■ 

the  Basils  mean  @K,  based  on  a  .-t-raple  of  sine  K  from  (say)  an  exponential 

distribution  with  mean  Q.  Far  hen  E(£L}8)  «  6,  all  Kx  and  11®  Var^te)  “0, 

&  K-ou  R 

from  which  the  above  condition  can  easily  be  shown  to  hold.] 

Then  we  have  the  following 

RESULT:  The  d.f.  H^T^)  *  P(Tk  <  t)  converges  pointwise  to  the  d.f. 

G(6)  *  P(8  <  t),  i.e.,V  fixed  t  >  0,  such  that  Hr  H„(t  )  «  G(t  ). 

°  *■ 

PROOF:  Fix  t  •  By  definition, 

-  — ■- G 


Tj  a 

w  “  /  /  g(6)fK(t|fi)dfi  dt 


c  o 


oo  MD 

/ s(i,)  /  fK(t|0)ett  dg,  assuming  the  ap*ropriate 


aeasurability  conditions  for  Fublni’s  Theorem  bold. 

Row  set  <^(6)  -  g(8)  J  fx(t|g)  dt,  sc  that  ^(tp)  - 


Since  v  d  and  X,  0^(6)  is  bounded  above  by  sup  g,  we  apply  toe 
dominated  convergence  theorem  and  get 

W  *  •  /‘!£  V9»*> 


•  f  Un  [«(«>  y*  fK(t|g)dt]  d& 

o  0 

/  A 

•  J  «(«)  Cil*  J  fK(t|g)dt]  da 
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SECTION  10. 0  COMPUTER  PROGRAMS  EEVELGiED 

10.1  DATA  ANALYSIS  PROGRAMS 

10.1.1  "KR" 

This  program  uses  data  in  pairs  consisting  of  S's  and  K^s  to  calculate 

the  moments  for  the  K-R  method  discussed  in  Sections  4.1  and  4.3,1  of  this 
report.  The  program  is  designed  such  that  if  either  the  second  or  the  fourth 
moment  is  negative,  it  terminates  and  prints  "method  not  applicable."  If  both 

of  these  moments  are  positive,  the  sample  82*  611(1  tae  estimated 

inverted  gamma  parameters,  o  and  are  calculated.  A  list  of  the  program 
follows. 


F?3p  MOP 


P«*/I5/S« 


HR 

IP  pero  y 
»  PP!NT"n*T i  TFT  NO.:"y 
30  PPINT 
40  LET  'll?* 

*p  let 
*0  LET  M!:P 
70  LFT  "Mr* 

*0  PE»0  F 

OP  PPINT"  ","MTPE","NC.  OF  FAILIVEF" 

100  PPIMT" . —  " 

UP  FOP  ,1  =  1  TO  y 

120  FEPO  T,N 

130  PPIMT  .l,T,N 

140  LET  Ml  :MI+  T 

150  LET  PP:H?+fTt2)*»V<|M-|> 

i*0  LET  itrK3+CTt3)*f»t2)/(*P-H)*<N-»-2)) 

170  LET  «M:M4+<Tti)*f  N?3  >/<<*♦!  )*<  N+?)*OM.J  >  ) 
|R0  NEXT  «! 

!P0  PPIMT’ . . . . - 

200  PPINT 
210  LET  M|:"l/K 
220  LET  M?rM2/K 
230  LET  *:M3/K 
240  LFT  M4:W4/F 
250  PPINT"*!:"  H 
2*0  PPINT"  MR:’  ? 

27(1  PPJNT’f  :”f 
?c(i  PFJNT" M»:" Ml 
LET  r? 

300  LET  VJr(1J-j4M|*Pp+p*(IH|  tj> 

SIP  LET  'Mrpj.j*MI*'lP+S*(|N|  t4) 

SEP  PPINT 

330  PPIPT"F»M»LF  MOMFNTF  APOUT  MEAN" 

34P  PPINT 

350  PPI  NT’FECONP:"'*? 

3*0  PPINT"  WIPrr"U3 
37P  PPI  NT"  FPNPTM:"  W4 
3PP  PPINT 

3o P  IF  U2«p  THEN  420 
400  IF  »4«0  THEN  420 
41 P  RO  TO  440 

420  PPINT"  METHOD  NOT  APPLICAPLF" 

430  RO  TO  530 

44P  |  ET  PI  sf  v?  t?)/(V2tJ) 

450  LET  El :F I  »,5 
4*0  LET  F2:«4/0>?t?) 

470  PP1NT"F4NPLE  PFTPIr’PI 

4P0  PPINTTPMPLF  FO  POCT  OF  PETAU’FI 

400  PPINT’FRMPLE  FFTA2="F? 

300  LET  K I :P  |*( (P2+3)  I?) 

510  LET  K2:4*(4*P?-3*P| )*< 2*P2-3*P 1 -*> 

520  PPI HT"FRMPLF  KAPPA r"KI /YT 
530  PPINT 

340  PPINT" . " 

330  PPINT 

3*0  PFI  NT"  TTI  MATEO  INV  RAM*A  PARAMETER?" 

570  PPINT 

500  LET  L:(2*MP-U1I  t?))/fM?-(MI  *2>> 

500  LFT  4:I»|*M -|  ) 

*00  PPI NT"ALPHA :"A 
*10  PPINT" LAMPPA:"L 
*20  PPINT 
*30  PPI NT*:::: 

*40  PPINT 
*50  RO  TO  10 
**P  EfT 


10? 


30.1.2  "BAYES" 

A  A 

This  program  first  calculates  a  andX.  using  equations  1.2. 2. 1.8  and 
4.2. 2. 1.9.  It  then  ufc.es  these  estimates  to  calculate  each  of  the  f(x.)'s 

2  1 

according  to  equation  4. 2. 2. 1.2.  The  X  cells  and  the  upper  class  limits 
are  read  into  the  program  and  it  computes  the  X*  value  by  equation  4. 2. 3. 1.1 
and  prints  out  a  XT  value  to  be  compared  according  to  the  method  described 
in  Section  4.2. 3*1.  Following  is  a  list  of  the  program  and  the  printed  re¬ 
sults  used  in  this  study. 


B*Y£5  *»3I  VEP.  19/f  J/60 

I*  P!H  Z<3#>,E(5*> 

Zf  PIN  P(2*),*<20),M<2*) 

3*  LET  Tt4J2* 

40  REAP  A 

*  PRINT  ’DATA  JET  NO. "A 

«  PRINT 

7*  LET  T0t0 

0*  FOR  1:1  TO  50 

3*  LET  Zd  >:* 

l«*  NEXT  I 

US  LET  %0 

12*  LET  St* 

IS*  LET  Yt * 

14*  REAP  X,F 

150  IE  Kt900  THE*  240 

168  IF  X»S  THEN  100 

IT*  LET  Z*tF 

1*0  00  TO  200 

10*  LET  Z(F)rF 

2*0  LET  Y*Y*F 

21*  LET  •fcPHF*  F 

22*  LET  Srf+(K»2)*K 

230  00  TO  14* 

240  LET  JfcH/y 
250  LET  StS/Y 
260  LET  VtS-(W2) 

27*  LET  Ar<l%T)/<V-M) 

200  LET  L*(  A*  N)/T 
20*  PRIRT  "ALPHA t”A 
3*0  PRINT  “  LANPPAi'L 
31*  PRINT 
32*  LET  PtT/(T+A) 

33*  LET  QiA/(T+A) 

340  LET  E*=«JtL> 

35*  LET  E*tE0«Y 
36*  FOR  Xtl  TO  5* 

370  LET  Ft(OtL)*(PtX> 

30*  LET  Pit | 

30*  LET  P2t  I 
400  LET  WtL+X-l 
41*  FOR  Itl  TO  X 
42*  LET  P|*P|*N 
43*  LET  P2tl*p? 

440  LET  WtW-l 
45*  NEXT  I 

46*  LET  E(X)tF0P|/P2 
47*  LET  E(X)tE(X)*Y 

40*  NEXT  X 
40*  PRINT 
50*  PRINT 
51*  REAP  Cl 

33*  PRINT  "*l*  W  CHI  S0Um  CELLS* "Cl 

ss*  print  L,"IT\ 

56*  FOR  Iii’to’ci-I . . 

37*  REAP  P<!> 

3*0  PRINT  I,P(!) 

30*  NEXT  I 
«*•  PRINT  Cl 
610  PRINT 
628  LET  e<D*Z0 
63*  LET  M(  I  )*E* 


W9 


M  POP  1:2  TO  Cl 
CM  LET  fi(I>:9 

CS»  let  H(I):» 

C7*  NEXT  1 

*99  TOP  J-|  70  P(l> 

CM  LET  0(|):9(|)+Z(J) 

WC  LET  K(I):P(I)+E(.» 

TIC  NEXT  J 

729  TOC  1:2  TO  Cl-I 

we  roc  j:Pt i-n+i  to  pen 

Tee  LET  0(1):0(1 )+;(.)) 

799  LET  N(2):H(I)+E(J) 

TM  NEXT  J 
779  NEXT  I 

we  let  si=e 
Tee  let  tisc 
see  roc  i =i  to  ci-i 
etc  let  st  z$  i+fl<n 

829  LET  TlrTl+HH) 
esc  cry  t  t 
849  LET  8(CI):Y-SI 
899  LET  H(CI ):Y*Tl 
8G9  LET  C2:9 

879  POINT  *(7»SE0VED’,*IXPCCTCP* 

809  POINT  * . - 

009  POO  1:1  TO  Cl 
N99  POINT  0(1) ,H( I) 

819  LET  C2:C2+((G(!)-II(I))  t2)/H(I> 

*9  NEXT  I 

ne  point 

849  POINT  TNI  S0UARE:9C2 
889  POINT 

SCO  POINT9::::::::::::::::::::::::::::::::: 

879  POINT 
809  POINT 
8*9  00  TO  *9 
1999  PATA  I 

1992  PATA  1,13,2,2,3,4,4,4.9,1,6,2,7,1 

1994  PATA  0,4,11.1,13,1,17,1,19,1,999,888 

I99C  PATA  4, I, 3,7, 2 

1999  PATA  I, II, 2.1 1, 3, (9, 4, 7, 9, 9, C, 4 

1*19  PATA  7,1,9,4,9,7.19,3.13,1,14,3 

1912  PATA  19, I, IC,3, 19, 1,22, I, 44, I 

1814  PATA  988,999 

I9IC  PATA  9,1,2.3,4,8,12,15,3 

1919  PATA  1,14,2,12,3, 13, 4, «, 9,3, 6,2 

1929  PATA  7,1,9,1,11,2,13,2,17,1,29, ',29,1 

1922  PATA  989,999 

1924  PATA  8,1,2, 3, 4,9, 4 

IC2C  PATA  1,19,2,11,3,8,4,8,9,1,6,2,7,3 

1928  PATA  9,1,9,1,19,1,12,1,999,999 

1939  PATA  A, |, 2,4.9 

1932  PATA  1,22,*, 7,3, 1,4,4, 9,2,8, I, 7,3 
1934  PATA  999,999 
1936  PATA  4, I, 2, 4, 8 

1939  PATA  I ,21 •C,V,3,9,4,2,9,9,6,4,7,3 
1949  PATA  9,2,9,1,19.1,999,999 
1942  PATA  9, I, 2,4,9. 7 

1944  PATA  1,19,2,19,3,9,4,2,9,3,6,3,7,3,8,4 
1948  PATA  19,4, II, I, 12, 3, 13, I, 19, I, 999, 999 
194*  PATA  9.1,3,7,19,9 
1999  PATA  I, 19,2,19,3,7,4,9,9,3, 6,3, 7,1 
1992  PATA  9,1,9,1,12,1,19,1,33,1,999,999 
1*94  PATA  4, 1, 3, 9 
8999  ENP 


DATA  SET  HO.  3 

ALPHA:  727.436 
LAWBPA:  .756319 


NO.  OF  CHI  S CHAPE  CELL?:  6 


CELL  NUHBEP  IPPEP  LIP1T 


1 

1 

2 

2 

3 

3 

4 

4 

5 

P 

6 

OPSEPVEP 

fryPECTEI 

14 

22.4572 

12 

6.63259 

13 

5.2156 

6 

4.19199 

7 

18.4516 

7 

10.051 

CHI  SQUAPE:  21. 

9927 

DATA  SET  HO.  4 

ALPHA:  3689.57 

LANPPAr  2.60495 

NO.  OF  CHI 

SOUAPE  CELLS:  4 

CELL  4UHBEP 

UPPEP  LIHIT 

1 

1 

2 

2 

3 

4 

4 

CBSEPVEP 

EFPECTEP 

11 

iS.ftll 

n 

9.I4P51 

12 

13.4987 

18 

12.4797 

CHI  SOUAPE:  1.31404 

SS3SS3?SSSSSSS*S**?S«SSSSSS«SSSSSS« 
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DATA  SET  NO.  5 

ALPHA:  75(14. 73 
LAWDA:  4.03(17 


NO.  or  CHI  SQUARE  CELLS:  4 


CELL  NUMBER  IPPER  LIMIT 


1 

1 

2 

2 

3 

4 

4 

CBSEPVED 

EXPECTED 

22 

15.9710 

7 

0.67011 

5 

10.3033 

6 

4.96470 

0<I  SQUARE: 

5.60010 

DATA  SET  NO 

.  6 

ALPHA:  2251 

.61 

LAMBDA:  1.73413 

NO.  OF  CHI 

SQUAPE  CELLS:  5 

CELL  NUMBER 

•PER  LIMIT 

1 

1 

2 

2 

3 

4 

4 

6 

5 

•SERVED 

EXPECTED 

21 

10.3603 

S 

0.79374 

10 

12.7930 

« 

7.3346 

7 

7.7*931 

OU  SQUARE*  I .30224 


:::::::::::::::::: 


U3 


!*TA  SET  NO.  7 


ALPHAS  |7«7.S9 
LAWPAs  •  »"7CC2 

NO.  OF  CM!  SQUAPE  CELLS:  9 


®Ll  NUPBEP  UPPEP  UNIT 


I 

1 

2 

3 

3 

7 

4 

If 

5 

CPSEPVEO 

EHPECTEP 

IS 

13. 2222 

19 

I4.P429 

II 

13.7991 

P 

€.23432 

€ 

4 .SNIPS 

OU  SQUARE:  9, €4293 

ssrsssrrsrssssr 

PATA  SET  NO.  p 

ALPHAS  770.SN4 

LANPPAs  .69229C 

NO.  or  CHI 

CELL  NUWPEP 

SQUARE  CELLS:  4 

UPPEP  LIHIT 

1 

1 

2 

3 

3 

9 

4 

QPSEPVEP 

DCPECTEO 

1C 

21.4934 

17 

If.fSS 

P 

S.II72S 

9 

(2.3434 

OM  SQUARE:  7.C7972 


3s:sztsssst;sssrss:sss;:ssss=sr 


Ilk 


DIMENSION  N(6),K(4),NCEL(6) 

DIMENSION  IPF (6,4,21 .TMETI200) 

01  PENS  I  ON  GPARKt 19I.CSI  (61  ,CS2(6I,Q(200I 

01  PENSION  7  PARK!  19I,NGTAI(20I.EG(201,ENI29)  .EPI20I ,0N(  191 ,0N(  191 
01  PENS  I  ON  INPFI6,4,2),INPF(6,4,2I,ICCLUM(4) 

01  PENS  I  ON  1(.U6),1C2I6).IC)(6> 

01  PENS  ION  XGAM(99l,XkE 18(99) ,XNRP(99I 
01  PENSION  VCAM(99,4| ,YP£I8(99,4I ,YNRM(99,4> 

DIMENSION  IUPF(6,4,2I,1C4(6) 

01  PENS  I  ON  EU ( 20 1 , DU ( 20 1 
IX*48828125 
5  FORPAT (6151 
7  FORPAT I  6 1 S ) 

IC  FORMA T ( 41 51 
IS  FORMAT ( 6F 10 . 3  I 

2C  FORMAT!  10F8.5/10F  8. 5/ IPF B. 5/ JOFfl . 5/ 10F 8 . 5 1 
25  FORPAt  (FIC.2.I  10.4F10.S) 

)C  FORPAT  (IHI/15H  INVERTEC  GAPPA.20X.7H  AL PHA* , 3X, F 10. 5 , 1  OX ,6H  LMDA* 
1.4X.I10  / 7H  L0GNRP,28X,4H  MU*  ,6X , F 1 0. 5 , IOX , 7H  S I GMA* , 3X , F l 0. 5/6H 

2WEIBLLL,?2X,7H  4LPHA* , 3 X . F 1 5.5 , 1  OX , 6H  BE T A* »F 10. 5/// ) 

35  FORMAT! 10X.3H  \* . I  4// 10X , 1H  K*,I4//> 

45  FORPAT ( //// ION  CELL  NO.  ,10H  CLASS  ,LOH  OBSERVED  ,10H  EXPECTED 
1.1CH  EXPECTED, 10M  EXPECTED  /10X.10H  MARK  , 10X.10H  INV  GAMPA, 

2  10H  LOG  N0RPL.10H  PtlBULL  //) 

5C  FORPAT(IIO,FIC.5.110,3F10.5) 

55  FORMA  T ( I  10, LOX, 1 10*3F10.5////20X, 14  H  CHI  SO  V ALUE S/ /20X . I  OH  INV  GA 
IMMA.10H  LOG  NCRPL.lOH  WEI8ULL  //OH  CCPPUTEC  CHI  SCUARE , 3F 1 0. 5/ 
220H  90  PERCENT  LEVEL  .3F10.5/2CH  95  PERCENT  LEVEL  .3F10.5) 

60  FQRPA Tl  //// 3CM  ***************************** f /// n 
TO  FORPAT  (///20X.25H  SUPPARV  OF  RESULTS  AFTER, 15, 8H  TRIALS/ 

120X.25H  LEVEL  t * . 90  LEVEL  2*. 95/  20X,  7H  ALPHA* , F 10 . 5 , 5K , 8H  LA 

2PB0A* , I  10/// ) 

75  FORPAT ( 20X, 7H  LFVEL^,|2///?0X,15H  INVERTED  GAMMA///) 

77  FORPAT  (10H  N*,6II0//I0H  K= / ) 

80  FORPAT  ( 71 10/1 

85  F0RPAT(///2GX,11H  LOG  NCRPAL///) 

87  FORPAT I ///20X.8H  PE  I  BULL // / ) 

88  FORPAT  (///20X.8H  UNIFORM///) 

9C  FORPAT ( ///20X , IQH  **•♦***♦*////) 

91  FORPAT  | 71 IC/8X,2H  L,6I10/8X,2H  *,6II0/«X,2H  U.6II0//I 

96  FORPAT  I///20X.12H  ALL  4  TESTS///) 

97  FORPAT  ( I5,F8.2,8X,FB.2,8X,F8.2) 

99  FORPAT  | 34H  TABIC  OF  SAMPLE  PERCENTAGE  POINTS/ 

12 JH  CF  THE  MARGINAL  FOR  X * , 1 5// 5X , 1 6H  INV  GAPM A  ,I6H  LOG  N.  RM 

2AL  ,RH  hE I  BULL // I 

142  FORMAT  I 10F8.2/IOF8.2/10F8.2/1OFB.2/10FB.2/I0FH.2/10FB.2/1OF8.2/ 
C10F8.2/9FB.2) 

RE  AC  (5,5)  (N(J),J*1,6) 


I 

i 


oc  in 

RE  AC  18,7)  (NCFl  (J).J=l,6t 

CC  1  i 

RF  AC  15,10  IK(J),J=l,4l 

ccso 

READ  (5,151  (CSl(J),J*lt6l 

CC4  i 

READ  (8, IS)  !CS?I  JI.JM.61 

cp4? 

ICC 

RE  AC  (5,75)  ALP.l AM.AMU.SIG.wAI  P.wBET 

CCA  t 

IF(ALP.E0.999.)  Uu  1 C  090 

CC'.'. 

CO  1C?  J=l,4 

0044 

RE  AC  ( 5 ,  1 42 i  ( AGAR  I  I 1,1  =  1,99) 

CC46 

RE  AC  (5,142)  ( XNRM  1 11,1*1,94) 

CC'.  7 

KEAC  (5,142)  ( XWC !B( 1 ) , 1*1 ,99) 

CC4M 

KO*(t  <  J  1 

CC<.  'i 

WRITE  (5,991  KC 

CC5P 

00  46  1=1,99 

OCsl 

WRITE  (6,971  1  ,XOA*(  I  )  ,  XNR*I  I  )  ,XWF  1  31  1  ) 

CCA? 

08 

CON  T  I NUE 

CCS  1 

WRITE  16,60 

CCS  <. 

no  l C  3  1  =  1,94 

rcss 

*G=XGAR( 1 ) 

r  356 

YGA7M  1  ,  J  )  =  XG 

CCS  7 

XN=XNMM( I ) 

ccsr 

YNRM ( 1 . J ) =XN 

ccs? 

<W=XWE I B ( 1 1 

COM! 

1 C  3 

YWEiet  l,J)  =  Xw 

CCfc  l 

IC2 

CONTINUE 

CCS.’ 

AUME*C. 

CON  i 

BUN  IF *2. *( <  J. /?.)**. 5 1*503. 

CC64 

WALP=(AMU/GAMMA< l.*l./WBET) )**WBET 

ccts 

WRITE  (6,301  AlP.lAM,AMU,StG,WAlP,wBET 

cot  6 

AMU l=AlOG(AKU**2/< ( AMu**2 *S I G**2 ) ** . 5 ) ) 

0067 

$  1 G 1 « ( ALOGI ( SIG*»2*AMU**2)/AMU**2) )»*.5 

CC68 

AMU=AMU1 

CCt  7 

S 1 G  =  S IG  l 

CC  7 P 

CO  /1C  1=1,6 

Cf  71 

CO  72C  J  = 1 , 4 

CO  7? 

no  73C  KK  =  1 , 2 

PC  7  7 

(PF 1 | , J,RK)=C 

CC  74 

INPF( I , J,KK  )  =C 

CC  74 

I WPF ( 1, J.KKl-0 

CC  77. 

IUPF'*  .  i-O 

nil 

73C 

CON  T 1  CUE 

CC7A 

72C 

CONTINUE 

CC  70 

7 1 C 

CONTINUE 

oceo 

1  TEN'C 

CCd  l 

no  7CC  LCPP=t,lCPP 

CC8? 

I  =6 

COR  ) 

NCE LL  =  NCEL (  I  ) 

CC84 

117 

NO«N( 1 1 

cops 

7N=NCFLl 

4 


I 


CC86 
ace  7 

CCf>8 
CCS9 
CC90 
OC<)t 
CC92 
0C93 
0C95 
0C95 
OC<56 
CC97 
0C98 
0099 
01CC 
0101 
010? 
01C3 
01C5 
01C5 
01C6 
010/ 
OlCfl 
01C9 
0110 
0111 
Oil? 
0113 
0119 
0115 
0116 
Oil? 
0118 
0119 
0120 
0121 
0122 
012? 
C125 
0125 
0126 
012/ 
0128 
0129 
0130 
0131 
0132 
013  3 


>*CELL*NCELL-1 

XNC*NC 

00  12C  J- 1  #  5 
KO*K ( J ) 
xkq*kc 

CO  12/  14*1.99 
XG*VGAN<  14,  J) 

XN*VNRPt IA, Jl 
XWVVlElBUA,  J1 
XGAfU IA1*XG 
XNRKI  IA)*XN 
XWE I B ( 1 A 1 «XW 
127  CONTINUE 

CO  119  JX«l.BCELl 
2J»JX 

ZMARK(JXl«ZJ/ZN 
Zf»«ZPARKl  JXI 
DO  123  L«1.99 
XL«L 

XL  »Xl.  / 1  CO  • 

IF  (ZN.GE.XU  GO  TC  123 

GM«XGAN(L-1  !♦(  l/H-XL*.ri)/.01>*lXC.AN(L)-XGAt‘tL-l> ) 
GO  TC  118 
123  CONTINUE 

118  GHARK1 JXI.GP 

119  CONTINUE 

DO  121  JX*l .NCELL 
U»GFARKUX) 

CALL  LOOKUP  (XWEIBtUrDhE) 

0M1 JXl«OhE 

121  CONTINUE 

DO  122  JX-l.NCEU 
U»GP4RK( JX| 

CALL  LOOKUP  ( XNftH , U f ONC ) 

ONI JXl-ONC 

122  CONTINUE 

00  750  JX«l,NCELL 
U«GPARK( JXt 

IF  ( U«  L T.bUN I F )  GC  TC  750 
DU( JX1-1. 

CO  TC  750 

75C  IF  (U.GT.AUNIF)  GO  TC  76U 

oui jxi«o. 

GO  TO  750 

760  DU1 JX 1*  I U-AUN1F  1  / 1 BUMF-AUNIF  1 
75C  CONTINUE 

EG1 1 1 •XNC»ZPARK 1 1 ) 

EGINCELL 1»XN0*1 1 • »2NARK INCELl 1 1 


0134 

0135 

0136 

on; 

0138 

cm 

0140 
0141 
014? 
0143 
0144 
0145 
0146 
0147 
014  8 
0146 
0150 
0151 
015? 
0153 
0154 
0155 
0156 
0157 
0158 
0156 
0160 
C  16  1 
016? 
0163 
0164 
0165 
0166 
0167 
016  8 
0160 
0170 
0171 
017? 
0173 
0174 
0175 
0176 
C  1  7  7 
0178 
0177 
0180 
C1P1 


EN(  1 ) *XNC*0N ( 1 1 
£N(NCELLI*XNQ*( 1  ,-DN INC ELL  1 1 
EM ( 1 I  «XNC*CH ( 1) 

LM(NCELU*XNC»<  l.-OW(MCELU) 

EU( 1 I  *XNC*OU ( 1 1 
EU(NCELLI*XNO*(  l  ,-OU 1 NCELL )  I 
IF  1 NCELL  >E  G . 2 1  GO  TC  126 
00  124  JX*2»NCELL 

EG(  JX)«XNC*(ZMA«Kl  JXi-ZMARKlJX-l 1 1 
EW( JX|*XNC*(0W( JXI-Oh( JX-lf I 
£NtJXI«XNO*ION(.JX)-ON(JX-l)  I 

tunxi«xNO*(ouuxi-oui  jx-i » t 

124  CONTINUE 

126  CO  125  11*1, NO 
01  1 1 1*0. 

125  CONfINUE 

00  130  11*1, NO 
CALL  AAN001 I  X, IV, ¥» 

I  X- I  V 

CALL  VRTGAN  ( ALP, LAN , V , THE  I A  I 
THE  flit  MTHETA 
00  150  JJ*l,KO 
CALL  AANDOMX.IV.ZI 
1 X- I Y 

S— THETA*ALCG(2) 

oim*Q<m«s 

15C  CONTINUE 

onn*oim/xKc 

1 3 C  CONTINUE 

00  170  1 T*1  .NCELL 

ngtali  i  n*c 

1  7 C  CONTINUE 

00  180  11*1, NC 
CO  ISO  IT*1, NCELL 
IF  (Ctm.LT.GNARKlini  GO  TC  200 
GO  TC  19C 

2CC  NGTALI  in*NGTALI  !T1»1 
GO  TC  180 
l«"  CONTINUE 

NGTAL (NCELL l-NGTAL (NCELL 1*1 
180  CONTINUE 
CSIG'O. 

CSN-C. 

CSM*C, 

csu»c. 

00  230  1T*1, NCELL 
F-NGTALI 1 Tl 

CSIG*CSIG*( (F-FG( ITI I««2I/1EGI IT) 1 


0182 

cm 

0185 
0185 
0186 
0187 
01  ua 
0189 
0190 
0191 
0192 
0193 
0195 
0195 
0196 
0197 
0198 
C 159 
0200 
0201 
02C2 
02C  3 
C2C5 
0706 
0206 
02C  7 
02C6 
02fc9 
0210 
C2U 
0212 
0213 
0215 
0216 
0216 
0217 
0218 
0219 
0220 
0271 
0722 
0221 
0225 
0225 
C226 
0227 
C228 
0729 


C$K*CSN5(  If-EM1TIS**2)/1EMI!T1) 
CSW*CSW+< IF-EWI l 1 1  l*»2)/ICWUTn 
23C  CONTINUE 

NGTAL(  191  *NGTALU9  J*NGTM.  1201 
EU( 191*£U{ 1 9 1 ♦£  U( 20 1 

231  00  2 32  I  T*l , 19 
f*NCMl(lTI 

csu«csu+< if-Eci ir i »**?»/ 1  eut in i 

232  CONTINUE 

IF  (CSIG.CT.CSK  I  II  GC  10  260 
(PF([,J,l>*lPF((,J,ll«r 
260  IF  (CSIG.GT. C$21  III  GC  <C  280 
IPFtl,J,2)*|PF(l,J,2>*t 
290  IF  (CCN.GT.CSll 1 11  GC  TO  300 
INPFI  1, j,  ll-INPFl I, J,l)*l 
3C0  IF  (CSN.CT .CS2(  I  I  I  GC  TO  320 
INPFI I,J,2I*INPF( I, J, 21*1 
12C  IF  (CSW.Gf.CSl  (II  I  GC  TO  350 
IWPF(  I.J.ll  -IWPFS  I, J, 11*1 
350  IF  1CSW.GT.CS2UM  GC  TC  765 
INPFI I,J,2)*IHPF( i,J,2l*l 
765  IF  (CSU.GT.25.9891  GC  TC  767 
IUPF(l,J,ll«tUPHi.J.l)M 
767  IF  (CSU.GT.28.869i  GC  TC  120 
IUPF{ I.J,21«1UPF( I, J, 21*1 
12C  CONTI NUF 

I  TEN*  I  TEN*  l 

IF  (ITEN.lT.50l  GO  TC  700 
WRITE  <6,70*  LCOP.ALP.lAP 
DO  37C  KK*l,2 
(.RITE  (6,751  KK 
WRITE  (6,771  (N(N»,N-l,6» 

00  380  J*  1 , 5 
DO  39C  1*1.6 
IC0LUP(I1*IPF(  l,Ji«r.| 

39C  CONTINUE 

WRITE  (6,801  MJ»,IICCluP(n,l*l,6> 
380  CONTINUE 

WR 1 TE (6,951 

WRITE  (6,771  (M*i,W«l,6» 

00  58C  7*1,5 
00  59C  (*1,6 
ICCLUNI  )*INPF(I,J,KKI 
59C  CONTINUE 

WRITE  (6,8(5)  K(J|,(|CCIUP<II.  1*1,61 
580  CONTINUE 
hR  I  T  E  1 6 , 8  7 ) 

WRITE  (6,771  ( N ( P I  * W* 1 ,6 1 


120 


0 


02  30 
0231 
02  32 
0233 
0234 
02  35 
0236 
02  3  7 
023e 
0239 
0240 
0241 
0242 
0243 
0244 
0245 
0246 
0247 
0248 
0249 
C25C 
0251 
0252 
0253 

0254 

0255 

0256 

0257 

0259 

0259 

0260 

0261 

0262 

0263 

0264 


00  580  J*  1 , 4 
00  590  1*1,6 
1C0LUW1  !  )=IWPF(!  ,  J.KK) 

59C  CON  f  1  NUE 

WRITE  (6,801  K(JI,(ICCCUW(I  )  ,  I  *  1 , 6 1* 

580  CONTINUE 

WRITE  (6,88) 

WRITE  ( 4 , 771  IN(W),M*1,6) 

00  880  J*l,4 
DO  840  1*1,6 
ICCLUN I 1 |*I UPf ( I, J.KK) 

090  CONTINUE 

WRITE  (6,80)  KiJ),IICCLUMI),l«l,6> 

880  CONTINUE 

WRITE  (6.96) 

WRITE  (6,771  ;N(M1,M*1,6) 

CO  680  J*l,4 
00  690  I*US 
!CI(I)*1PP(!,J,KK) 

IC2( I >* INPf ( ( , J, KK  I 
IC  3  ( I  I  •  I  WPF  (I, J.KK) 

IC4( I J*»UPEI T.J.KKI 
690  CONTINUE 

WRITE  (6,41)  K(JI,IICl(L),L*l,6),(IC2(L),L*l,6),(IC3(L),L*lt61 
C.i IC4(U,L*1.6) 

680  CONTINUE 

UR  I TE ( 6, 40 1 

370  CONTINUE 

WRITE  (6.642)  IX 
692  FORM!  (10X.4H  IX*,  120) 

WRIIE  (6,60) 

I TEN'C 

7C0  CONTINUE 
GO  TC  100 
994  STOP 
ENO 


121 


COOl 

000? 

ccc? 

0CC4 
0005 
0QC6 
0CC7 
00' 8 
OCCP 
00 10 
ecu 
0012 


SU0RGUT INE  LOCKUP  (ARRAY, U.P) 

DIMENSION  ARRAY<<J<11 
DO  ICC  L-1,<J<J 

I f  ( L .GE  AR RA  Y ( l ) I  GO  TO  100 

Pl«l-1 

Pl»Pl/lCO. 

P*Pl*»Ot*(U-ARRAY<L”lH/(  ARRAY  ( L  1  -  array  ;  L  - 1 1 1 
CO  tc  200 
ICC  CONTINUE 
2CC  CONTINUE 
return 
ENC 


122 


I 

i 

| 

V 

V 

| 


I 


i  CCCl 

f  ccc? 

COC  1 

l  ccc*. 

i  CCCS 

l:  OCCb 

i  OOC  J 

ccc» 


SU6RCLT Iftf  «»NCO  IIX,lY,YFll 
ItllT)  5.6,f 
6  VfflY 

VFL«m*.4A%66|1l'-<) 

RE U.RN 
f  NC 


CCCl 

CC02 

0003 

COCA 

CCC5 
CCCfc 
0CC7 
OOCE 
COCO 
COlO 
OOlt 
0012 
CC 1  3 
CCtA 
CCI5 
OClfr 
CC  l  7 
CC  1 H 
0019 
0020 
CC2i 
0022 
0021 
0C2A 
CC25 
0026 
0027 
0026 
CC29 
CCJO 


SUBROUTINE  VKTGAPI ALP.LAN.P, THE  I* I 
ALAP-LAM 

AMU«AUP/| ALAH-l.  ) 

50  XIM./ANU 
N*LAP-1 
ICCLNT*C 
ICO  V*  l  • 

S*l. 

00  ISO  K*1.N 
FK»K 

¥«V*(ALP*X1 >/FK 
S«S*V 

15C  CONTINUE 

X3*X2*XI S*EXP<_4Lf,*X1 *  l*“F* 1/1  *I-P*V*EXPI-ALP*X1 )  1 

FARS*ABSI X  3 1 

IF  (FABS.LT. 0. 000001  I  GO  TO  200 
X1*X2 

ICCLAT*1CCUNT*1 
IF  ( ICCUNT.EQ.20I  GC  TO  250 
GO  TC  ICO 
200  T *  1 .2X2 

IF  IT.GT.O.)  C.C  TO  300 
25C  t*R  I  TE  1 6 >  2A0 1  P 

2A0  FORPAT  I//12H  HUNG  UP  ON  .F10.5//) 

260  P«RAND( X I 
GO  TC  50 
3CC  THE  TA*  T 
RETURN 
END 


* 


124 


oooi 

DIMENSION  V«T:.15ll  ,WF  IM511  t  2LUG13  1 1 

OOO  2 

Of  MEN  SION  IRANIM  511  ,X  VR  T5  151 »  ,X  HE I  8  {51 1  ,<  ZLOGf  5 1 1 

00".? 

DIMENSION  1J(  201  ,10.(2  01  »MG(20I»HM(20I,-(L(20I 

0004 

DIMENSION  C  UMI  3  (  51 1 

00"  0 

I X*48S  281 2  5 

000ft 

ICO 

7  F  AD  I  5,251  A  LP  ,A  I.AM,  AMU,  SIl. ,  NALP  ,  M  BEf  ,7  1  ME 

ooo  r 

25 

FORMAT  (  7F  10.  51 

000  8 

IF  IA1  P.EQ.999.  1  3C  TO  999 

000*» 

Ha  1TF  I  6,941 

0010 

54 

0011 

MALP*<  AMU/GAMMA  1 1.41.  /HUE  T!  1  ••MBET 

0012 

MR  ITE  (6,30  A  L  P,  A  LAM,  AMU  ,  SIG  ,MALP  ,d8E  T 

001  3 

30 

FORMAT  (  10X,2F10.5//10X,2F1C.5//UX,2F13. 57/71 

0014 

AM  L  1*5  LOU  I  A  MU*  *2  /(  I  AHU»A?  ♦  SI  G**2 1  M*.5 1  1 

0015 

SIG1*(  AL3f,<(  SIG**24AMU**2)  7AMU»»21  1  **.  » 

001ft 

AMI*AMU1 

0017 

sig*sig  1 

001ft 

00  in  1*1,  R 

0019 

VR  TGI  1  1*0. 

0020 

MF  ini  I  1*0. 

0021 

* 

110 

2L0G1  I  1*0. 

0022 

2  G=0. 

0023 

1  M*0. 

0024 

2L«0. 

c 

MAlGINAi  FOR  INV  GAMMA  PPM™ 

0025 

l 

P * T IMr  71  T1ME4ALPI 

0026 

j-alp/(  tim:»alpi 

0027 

VR  TGI  1  1*0*«ALAM 

0028 

XTAL* VA  TGI  11 

0029 

C  UMIGi  ll-XTAL 

0030 

on  120  J*1  ,7 

0011 

F»(Q*«UAM»*(P**J! 

00)2 

P  1*1. 

0033 

P  2*  l. 

0034 

XJ*J 

0035 

W*ALAM  *  XJ-  1. 

0036 

no  no  i*i, j 

003  7 

P  1*0  l*M 

00  3  5 

XI*I 

0039 

P  2*XI *P  2 

0040 

W»krl. 

0041 

no 

rONTINUT 

0042 

VR  TG(  J*l>«f  *P1/P2 

0043 

XTAL*XTAL4VPTG(  >1) 

0044 

CUM1GI  J»ll*XTAL 

0045 

120 

CONTINUE 

126 


004d 

r 

7G*l.-XT4l 

K 

c 

r 

P4#GIN4l  60#  MFIrtUll  PRIOR 

004  7 

L 

00  200  1*1,10000 

0048 

CALL  AANOOIt  X,:  Y,#l 

0044 

IX*  1 Y 

005J 

THf  t»*l  MLP*AUC(1./PII  ••Cl./METI 

0041 

H*0 

0052 

CUH"k«0» 

0053 

2  5" 

CALL  A4ND0M  X,l  Y,M 

0054 

1 X*  l v 

0055 

T»-Tltr  T4  *A  1 03 1 r  1 

005?. 

elec  «:lo:x*t 

0057 

irCLor.K.3T.Tf  Hf  1  SO  TO  300 

0058 

IHK.L  T.8I  30  TO  2  01 

0054 

7  8*7  8*1. 

00  50 

GO  TO  200 

00c.  1 

201 

.  *X>  1 

0062 

SO  TO  2  50 

0063 

3on 

ME  l8IKM)*«l8(K»litl. 

0064 

200 

CONTIMUf 

0065 

JO  110  1*1,8 

0066 

310 

Mb  (81  1  1**18(11/10000. 

006  7 

r 

iM*7M/lOOOO. 

l 

c 

c 

MA#r.lN4L  60#  LOG  NO#  UAL  PRIOR 

0066 

L 

00  400  1*1,10000 

0060 

A*C, 

0070 

00  450  <1*1,12 

0071 

CALL  A*  XDO( !  X , I  Y, Yt 

0072 

I  X*1  Y 

0073 

450 

4*4*Y 

0074 

T«ET4.FXP(  14-6.  01  *SI344PUI 

0075 

x»0 

00  7  o 

CLPCX-O. 

0077 

46  0 

CAIL  AARJOU  X,l  Y,F| 

0078 

I  X*  (  V 

0070 

T*-TH?  T4  *4  10GI  6  1 

.7080 

clocx*:io:  x*t 

0081 

IF(CLJCX.3T.TIPEI  GO  TO  480 

0082 

If  5X.LT.81  30  TO  470 

008? 

7L*7l*l. 

0084 

GO  TO  400 

0085 

4  70 

K*X4l 

0066 

SO  TO  460 

0087 

480 

7LCG(  #♦  1 1  *71031  KM»*1. 

OCttf 
006  9 
0090 
009? 
009? 
009  3 
009* 
0095 
OOVo 
009  r 
0098 
0099 
0101- 
0101 
0102 
0103 
0109 
0105 
0106 
0107 
0106 
010“ 
0110 
Olll 
0112 
0113 
0119 
0115 
0116 
OUT 
0118 
0119 
0120 
0121 
0122 
0123 
0129 
0125 
0126 
0127 
0126 
0129 
0130 
0131 

on? 

013  3 
0139 
01  35 


9C0  CONTINUE 

00  910  !«? ,u 

*10  2 LOG!  1  1  *21031 !  I  /I  0000. 

U’ZL/ 10000, 

01  500  1.1,6 
bE*0(5,5|  NO 
5  FORMAT!  151 

•*  ITF  I  6,961  NU 
96  FORMAT  1//10X.3H  N*,f5) 

* t*0*  5,101  i:i.CSl,CS2 

io  haunt T(  r  io,2f  io.  6i 
jc  t*  i:  i- 1 
00  too  IM.JC1 

9EAni5,?|  101  I  ; 

101  1  1*101  1  »*1 
(Or  CONTINUE 
XND«NO 

00  TOO  1*1,8 

xvp  *r.i  ii*vPT3in*xNn 
*"H8i  n*wtiftin*xNn 

X210GI  I  )*2lrx?(  1 1  *XNU 
TOO  CONTINUE 

x2g*/s*xnu 

XI  6*2  W*XNO 
Wl  *21  *XNO 
I8Cl*0 
IP*1«0 
IPt  1*0 
IPG2-0 
IP  62-0 
lpl  2*1 

nn  90i  mup. i,i0oo 
DCS  990  1*1,8 
950  IRANI!  I  1-0 

u*o 

UO  610  I  *1  ,ic  1 
IG1  |  »*o 
MG1  1  1*0. 

HMIU-O. 

HtllJ-O. 

6JC  CONTINUE 

DO  511  11*1, NO 
C*U  >*N30llX,|  Y,2» 

I  X*  I Y 

DO  520  J-1,6 

IT «l T.C UMT3I J) |  50  To  530 
*20  CONTINUE  0 

18  *!»♦  l 
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go  ro  510 

530  |R«N0(  JI*|FAN0U1»1 
510  CONTINUE 

IF  (LOOP.Nf.il  GO  TO  01  0 
WR  ITFt  5.991 

99  FORMAT  (///18X.2H  K.5X.9H  FXPLCTED.6X.1H  EX  PECT  EO,  6K,  9H  EXPECTED, 
16X.9H  08£P^O/?5X,10H  I NV  GAMMA, >X  ,10 i  MEIOJLL  ,  5X,  10H  LUG  NOftNL 
2// I 

00  800  1*1,8 
lJ*t-l 

WRITF  (6,981  I  J.XVRTGd  I  ,XWLI6(II  ,X2LUG(  U.IRANOU  I 
98  F3PMA  T(  10X  ,1  1  0, 5X ,F t  0.  5  ,5 X ,F 1 0.  5  ,» X  ,Fl ) .5  ,5«  ,  1 13  1 
800  CONTtNUT 

WRITE  (  6,9T|  *16  ,X2W«K2L«I  R 

9T  FORMAT!  10X.8H  OVER  T  , TX  ,Fl 0.  5  ,5X  ,F  10.  i  ,5X  ,F1>  .5,5X,  1 101 
910  !l)l«n(  1) 


0182 

0183 

0184 

0188 

0188 

0187 

0188 

01u9 

0190 

0191 
019? 
0193 
0194 
0185 
0196 
0197 
0198 
0199 
0200 
0201 
0202 
0203 
0204 
020  5 
0206 
0207 

0208 

0209 

0213 

0211 

0212 


00  680  i»i  .in 

G»  IGt 1 t 

C  SG»C  SG*(  (  G-HGI  III **2 1 /HG II ) 

CSk»CSM«((G»4Hmi  »*2)  /Mali) 

cst*c  si»iirHimt«n/HUii 

680  CONTINUE 

If  (t3nP.NE.il  GU  TO  920 

tft  ITFI  6,93)C  SG  ,C  $W,C  SL  ,CSl  ,CSl  ,CSt  jCS2  ,CS2  ,CS2 
93  FORMAT!  ///10X.18H  CNi  SQUARE  VALUES// 1 )  <  .  3F 1?  .5/ 13X,  3F 10 . 5/ 10X.3F  1 
C  0.5//) 
i*t  IT£  I  6,911 

95  FORMAT  1///19X,10N  **•***•*»////) 

920  IF  (ESG.GT.CSII  GU  TO  930 
1PG1-1PG1M 

930  IF  (CSG.GT.C  12)  GU  TO  940 
IPG2-IPG2*! 

940  IF  ICSK.GT.CSll  GO  TO  950 
!P*l"!PWl»l 

550  IF  (CSM.ST.C  S21  GO  TO  960 
IP  K2*IP  H2*  I 

560  IF  (CSL.GT.CS1I  GO  TO  976 
IPL  1*IPL  1*  1 

970  IF  ICSL.GT.C  521  GO  TO  900 
IP!  2MPL2V1 
900  CONTINUE 


<«lTf  16,921  1  PG1  ,1  PG2  .1  PW1  ,1  PK2 ,1  PL1.I  PL2 

92  format  i///iox,0h  sonmby7/7i ox  7i  vx , 2 nn 


C10X.10H  KEIBULL  ,5X,2I  5//1  0X,1  OH  LOG  N09NL.5X  ,2  15// 1 


kRITE  16,941 
500  CONTINUE 


GO  TO  100 


999  STOP 


ENP 
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000!  SUBROUTINE  RilOO  1 1  X,l  V.YFLJ 

0002  ir«ix*6553S 

OOOJ  IFIIY!  5,6,6 

0004  5  ir«rv+2I4748J64M 

0005  6  VFl*  1 V 

0006  Yfl*VM.6.4fr566U£-<» 

0007  RETURN 

00"8  fcNO 
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ences  used  for  this  study,  whereas,  Section  11.2  is  a  general  bibliography 
containing  sources  of  material  on  Bayes  reliability,  in  general,  and  is 
meant  only  for  the  readers  use  and  convenience . 
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